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We develop a kinetic equation description of Coulomb drag between ballistic one-dimensional 
electron systems, which enables us to demonstrate that equilibration processes between right- and 
left-moving electrons are crucially important for establishing dc drag. In one-dimensional geome- 
try, this type of equilibration requires either backscattering near the Fermi level or scattering with 
small momentum transfer near the bottom of the electron spectrum. Importantly, pairwise forward 
scattering in the vicinity of the Fermi surface alone is not sufficient to produce a nonzero dc drag 
resistivity po, in contrast to a number of works that have studied Coulomb drag due to this mech- 
anism of scattering before. We show that slow equilibration between two subsystems of electrons of 
opposite chirality, "bottlenecked" by inelastic collisions involving cold electrons near the bottom of 
the conduction band, leads to a strong suppression of Coulomb drag, which results in an activation 
dependence of pr> on temperature — instead of the conventional power law. 

PACS numbers: 71.10.Pm, 73.21.Hb 



I. INTRODUCTION 

A remarkable property of a system of two conductors 
placed in proximity to each other is the occurrence of the 
phenomenon of Coulomb drag. This consists of induc- 
ing an electric field or current in one of the conductors 
by sending a current through the other — with the fric- 
tion force being due to electron-electron interactions — in 
the absence of transfer of electrons between the two sub- 
systems. As such, Coulomb drag is a sensitive probe of 
electron-electron correlations and, specifically, of inelas- 
tic electron-electron scattering. 

The key quantity describing friction is the drag resis- 
tivity pd defined for two homogeneous conductors paral- 
lel to each other as 

Pd = -E 2 /ji , (1.1) 

where j\ is the electric current density in ( "active" ) con- 
ductor 1 and E2 is the electric field generated in ("pas- 
sive" ) conductor 2 to compensate for the friction force un- 
der the condition that no current flows in the passive con- 
ductor. Since its prediction^ a third of a century ago, 
for two-dimensional geometry of two parallel conducting 
sheets, Coulomb drag has been extensively studied exper- 
imentally in double-layer semiconductor structures j 3 - - — 
also in a transverse magnetic field J£r— Recent exper- 
imental work has addressed a similar phenomenon in 
double-layer grapheneJ^— In one-dimensional geome- 
try, a number of experiments have explored Coulomb 
drag between quantum wiresi^Sr— Drag experiments 
have also been done on electron systems of other geome- 
try: between two- and essentially three-dimensional elec- 
tron systems 2 ^ or between quantum-point contacts^ 

A great deal of understanding of the mechanism of 
Coulomb drag has been achieved by calculating the fric- 
tion force perturbatively, at second order, in the dynam- 
ically screened interaction Vi2(oj,q) between two two- 
dimensional electron systems ("orthodox theory")! 2 ^— 
Within this framework, Coulomb drag is represented as 



rectification of nonequilibrium current fluctuations in- 
duced in the passive layer and, consequently, the linear- 
response resistivity pd is related to dynamical correla- 
tions in thermal fluctuations of the electron densities in 
different layers at equilibrium. Equivalently, pd within 
the orthodox theory is proportional to the rate of mo- 
mentum transfer between the layers at order V 2 2 - One 
important result of the orthodox theory is that pr> at or- 
der scales with temperature T in the limit of small 
T as T 2 (for ballistic electron systems, or as T 2 lnT in 
the diffusive limit). The power-law vanishing of pu as T 
decreases is associated with the constraints on the phase 
space available for inelastic electron-electron scattering. 

To the best of our knowledge, in all works where 
Coulomb drag in two-dimensional systems was studied 
within the framework of the orthodox theory, pd was 
derived under the tacit assumption that the intralayer 
relaxation processes (determined, e.g., by disorder) are 
faster than the processes of momentum transfer between 
the layers. Within the kinetic equation approach, which 
we employ in this paper, this condition implies that an 
iterative solution 29 of the kinetic equation (equivalent, 
in the diagrammatic language, to the evaluation of the 
Aslamazov-Larkin-type diagram a 31 i 32 ) is justified. A del- 
icate point here is that the resulting drag resistivity in 
the presence of disorder does not necessarily depend on 
the strength of disorder, which might seem to imply 
that the thus obtained pu describes the clean limit as 
well. However, in the absence of relaxation processes in- 
duced by disorder or inelastic intralayer interactions, the 
nonequilibrium part of the electron distribution function 
is governed by interactions between the layers, so that 
the lowest-order expansion in the interlayer collision in- 
tegral, assumed in the orthodox theory, is generally not 
sufficient. Therefore, the orthodox theory should not be 
expected to be generically valid in the clean limit — even 
for an arbitrarily weak interaction between the conduc- 
tors. One particular example that demonstrates a dra- 
matic departure from the orthodox theory in the clean 
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case is Coulomb drag between ballistic quantum wires, 
addressed in this paper. 

The orthodox theory also explicitly points to the im- 
portant role of electron-hole asymmetry in a degenerate 
Fermi gas, in the absence of which the electron and hole 
contributions to pd at order V 2 2 cancel each other. The 
cancellation^ 3 , has the consequence that, in the case of 
particle-hole asymmetry produced by a hnite curvature 
1 /to of the electron dispersion relation, where to is the 
electron mass, po is small in the parameter (T/ep) 2 with 
ep being the Fermi energy. In the diffusive case, po for 
sufficiently small transferred momenta can be directly re- 
lated at order V 2 2 to the dependence of the local con- 
ductivity on the local electron densit y 34 ' 35 absent in the 
particle-hole symmetric case. Apart from a nonzero cur- 
vature I/to, particle-hole asymmetry can also result from 
the energy dependence of the electron density of states 
in the vicinity of cf ■ The latter contribution to Coulomb 
drag is important in the presence of a transverse mag- 
netic fiel d 35 ' 36 because of the modification of the density 
of states by Landau quantization. It is also important 
in two-dimensional electron systems with a linear disper- 
sion relation (1/m = 0), in which the density of states 
varies linearly with energy; in particular, in graphene. 
Particle-hole symmetry in two dimensions is realized in 
graphene at the charge-neutrality point, where the or- 
thodox theory gives zero po- Therefore, possible devi- 
ations from the orthodox theory in the vicinity of this 
point in graphene become particularly important. 38 The 
prevailing notion that Coulomb drag is entirely due to 
particle-hole asymmetry is justified in the case of a dis- 
ordered two-dimensional electron system only at order 
V 2 2 ■ Beyond the golden-rule level, already at third order 
in V12 , rectification of interaction-induced current fluctu- 
ations in a diffusive double-layer system yields nonzero 
Pd even in a particle-hole symmetric system.— 

In one-dimensional geometry, the connection between 
Coulomb drag and particle-hole asymmetry is subtler. 
Processes of electron scattering due to interwire interac- 
tion separate into two classes: backscattering, in which 
an electron changes its chirality, and forward scattering, 
in which it does not. In a ballistic system with no disor- 
der, the contribution of interwire forward-scattering pro- 
cesses to Coulomb drag vanishes if the electron dispersion 
relation is linearized (Luttinger-liquid model 4 ^); however, 
the contribution of interwire backscattering processes to 
Pd remains nonzero even in the particle-hole symmet- 
ric limit Much of the prior work on Coulomb drag 
between quantum wires has therefore focused on the 
backscattering processes. At the golden-rule level, po 
induced by backscattering between identical wires is lin- 
early proportiona l 50 ' 51 to T with 



where /?b is the dimensionless coupling constant describ- 
ing interwire backscattering at the Fermi level with mo- 
mentum transfer 2fcp, vf is the Fermi velocity, and 



2ir/e 2 is the resistance quantum (here and below h = 
1). At higher orders in the strength of interaction, 
both intra- and interwire, a power-law renormalization 
of the backscattering amplitude develops as T — or the 
drive current in the active wire in the nonlinear re- 
sponse regime — decreases 4 ^ -45 ' 47 ' 49 ' 52 (a similar renor- 
malization of pd in the strongly-interacting limit of a 
"spin-incoherent" Luttinger liquid has been discussed in 
Ref.|H). 

Below a characteristic energy scale (at which the renor- 
malized amplitude <?i is of the order of unity), electrons in 
two wires form a "zig-zag ordered" charge-density wave 
and the power-law behavior crosses over into an exponen- 
tial growth of pd with lowering T or, in finite-size sys- 
tems at sufficiently low T, into an exponential growth of 
the drag resistance with increasing system size^ 2 - -44 ' 47 In 
the limit T — > 0, the resistivity po (defined as the linear 
resistance per unit length under the condition that the 
size of the system is made infinite before any other limit 
is taken, in particular, that of zero T) is infinitely l arg e 
( "absolute current drag" in the terminology of Ref. |42|) . 
By contrast, the linear drag resistance between finite-size 
wires vanishes to zero as T 2 in the limit of low T, inde- 
pendently of the strength of intrawire interaction and on 
whether the wires are long enough to form the zig-zag 
order or notj 44 ' 53 In the former case, however, there ex- 
ists a parametrically wide range of T in which the drive 
and drag currents are almost equal to each other (almost 
absolute current drag) up to an exponentially small un- 
balance due to transport of solitons in the charge-density 
wave. In the nonlinear response regime, the drag re- 
sistance is finite at zero T and shows oscillation o 42 ' 49 
as the drive current varies, related to the interference 
of plasmon waves reflected from the boundaries of the 
wires, which are suppressed^ 9 , as T is increased. If the 
length of the region in which interwire backscattering 
occurs is much smaller than the total length of the wires, 
Coulomb drag in the limit of low T and small bias volt- 
ages can be described in terms of backscattering at effec- 
tively pointlike contact] 41 ' 45 ' 46 ' 54 In this model, the drag 
conductance induced by electron-electron backscattering 
is expressibl o 41 ' 45 ' 55 in a particularly simple form in terms 
of the conductance of a single wire with a single static 
backscattering impurity. 

One of the conclusions that follow from the above re- 
sults for backscattering-induced Coulomb drag is that — 
even if the bare (before the renormalization) backscat- 
tering amplitude is small — the drag effects can be strong 
in the infrared limit, which for the linear drag resistivity 
Pd means sufficiently low T. However, the backscatter- 
ing amplitude falls off with increasing distance a between 
the wires as exp(— 2kpa), where Uf is the Fermi wave- 
length (assuming the electron densities in the wires to be 
equal to each other). As a result, for kpa ^> I the drag 
effects that originate from backscattering are strongly 
suppressed unless T is exponentially small and electrons 
are zig-zag ordered. Moreover, the effect of backscatter- 
ing is also suppressed in the case of nonidentical wires 4^ 
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An alternative mechanism of drag is interwire scattering 
with small momentum transferor— Despite relying on 
electron- hole asymmetry (e.g., a nonlinear dispersion re- 
lation for electrons) or the presence of disorder, 59 this 
mechanism of drag is expected to be more effective than 
backscattering if quantum wires are sufficiently separated 
from each other, and is further favored by an imbalance 
in the electron densities. 

Against this background, it is desirable to explore the 
possibility of Coulomb drag due to interwire forward scat- 
tering in the absence of any backward scattering. An 
important advance in this direction was made in Ref. [5(j| 
which extended the orthodox theory 2 ^— for electron sys- 
tems with a parabolic dispersion relation to one dimen- 
sion (see Ref. for a review of one-dimensional non- 
Luttinger liquid models) — under the assumption that 
electrons are ballistic (no disorder) and can only ex- 
change momenta much smaller than kp. Specifically, the 
result obtained in Ref. [56| for po induced by forward scat- 
tering in the vicinity of the Fermi level for T -C vf/o. 
reads 

p D ~ — f3 2 k F ( — J , (1.3) 

where /3f is the dimensionless coupling constant describ- 
ing interwire forward scattering. The meaning of the 
question mark will become clear in the next paragraph. 
At higher T, in the interval vf/cl <C T <C cf (which ex- 
ists for kpa^> 1), pd shows a plateau in the dependence 
on T according to Ref. Other wor k 57 ' 58 has reached 
similar conclusions by employing a bosonic description of 
the one-dimensional electron liquid with a nonlinear dis- 
persion relation for electrons; in particular, reproduced 
Eq. (jl.3|) . According to Ref. H3, however, the T 2 scal- 
ing of pd for the case of identical wires is only valid for 
/3fep -C T -C vpja (provided that vf/cl is larger than 
PftF and such an interval of T exists). The energy scale 
ftiCF describes splitting between symmetric and antisym- 
metric plasmon modes in the double-wire system. In the 
low-T limit, for T <C f3{£F, the drag resistivity between 
identical wires is predicte d 57 i 61 to vanish with decreasing 
T as p D - /3f 1 (2Tr/e 2 )k F (T/e F ) 5 . Importantly, all the 
prior work 5 - ^ 58 ' 62 on Coulomb drag due to interactions 
with small momentum transfer obtained a nonzero con- 
tribution to pd (unless T = 0) from forward scattering 
in the vicinity of the Fermi level. 

One of the purposes of this paper is to demonstrate 
that in fact 

Pd(T) = (1.4) 

in the absence of scattering that changes the chirality 
of electrons. That is, forward scattering itself cannot 
establish nonzero dc drag. As will be shown, the ortho- 
dox theory^Sr— with the use of which (or of its equiv- 
alent in the bosonized formulation of the problem) the 
nonzero result was obtained in the earlier works ^ r 58 i 62 
fails entirely — at the conceptual level — in one dimension. 



The basic question behind this problem is under what 
conditions the second-order expansion 2 ^— of pd in pow- 
ers of V12 , which constitutes the essence of the orthodox 
theory, is justified. Clearly, this is correct if drag is suffi- 
ciently weak. The question is about how weak. The an- 
swer, as we will demonstrate below, contains a delicate 
but crucially important point which does not appear to 
have been generally appreciated in the literature. 

The Kubo-type formula for po derived within the or- 
thodox theory in one dimension reads: 

Pd = ttt~ — m fw~ f p-1 2 \Vi 2 (u},q)\ 2 
2e i n\n2± J 2ir J 2tt 

Imn 1 (w,g)Imn 2 (w,g) 

x ^ , 1.5 

sinh 2 (o;/2T) 

where IT^w,*?) and n\ t i are the (retarded) polarization 
operators and the electron densities, respectively, in wires 
1 and 2. The product of the imaginary parts of ili^w, q) 
results from the application of the fluctuation-dissipation 
theorem to the equilibrium dynamical structure factors 
for density fluctuations S%(oj, q) and S2(— oj, —q). The le- 
gitimacy of the use of the lowest (second) order in V\2 for 
Pd is based on the assumption that the density fluctua- 
tions in the active conductor (wire 1) are equilibrium in 
the frame moving with the drift velocity Vd — — j\jen\ 
(throughout the paper the electron charge is defined as 
— e, i.e., e > 0). Indeed, as demonstrated in Ref. [H, 
the expansion of S\(uj — qvd,q) to first order in ji in 
the cross-correlation function of the electric forces in two 
conductors gives^ the linear-response dc electric field 
E% = —pr>ji in wire 2 (for j2 maintained at zero) with 
p D from Eq. (JT5J. 

However, the assumption about the density fluctuation 
being equilibrium in the moving frame ("drift ansatz") 
is by no means innocent: actually, in one dimension, it 
strongly limits the applicability of Eq. (jl.5[) . If one as- 
sumes, for dcfinitcncss, that the wires are identical (with 
the electron density n) and represents po as 

p D = m/e 2 nr D (1.6) 

by introducing the "drag rate" 1/td which describes 
momentum exchange between two conductors in the dc 
limit, the drift ansatz is only legitimate — as will be shown 
below — if 1/td is much smaller than the thermalization 
rateO In quantum wires, thermalization means not only 
energy relaxation within the same chirality branch but 
also "right-left" relaxation. The latter, however, can only 
occur if backscattering is allowed, so that in the model of 
Refs. IMIlll in which there is no backscattering "by 
construction" , the use of Eq. (|1.5[) is not legitimate. At 
this point, one might think that a deviation of the ex- 
act shape of the distribution function of electrons in the 
active wire from equilibrium in the moving frame does 
not change the result (|1.3p qualitatively, i.e., only the 
numerical coefficient in Eq. (jl.3|) depends on the shape 
but remains of the order of unity. This is, however, not 
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the case; on the contrary, the shape is exactly such that 
in the absence of backscattering p D is zero [Eq. (|1.4p ]. 

Naively, one might expect that the results of Refs. [56l — 
for drag induced by forward scattering are valid for 
sufficiently long ballistic wires, namely for wires whose 
length is much larger than the characteristic scale of the 
right-left relaxation. We emphasize, however, that the 
equilibrium state in the moving frame in the double- wire 
system cannot be reached by increasing the length of the 
wires beyond this scale (as would be the cas o 76 ^ 80 in a 
single wire). Nonzero friction prevents this from hap- 
pening by constantly exciting electron-hole pairs — even 
in the limit of an infinitely large system size. 

Our main results can be described as follows. We 
demonstrate that pr> in one-dimensional geometry van- 
ishes in the case of electrons not changing their chiral- 
ity in scattering processes. A key consequence of this 
is that Coulomb drag between quantum wires is neces- 
sarily suppressed compared to the conventional theories 
[epitomized by Eq. (|1.5|) ] if the right-left equilibration is 
not fast enough. In the case of energy relaxation being 
mainly due to processes with momentum transfer much 
smaller than kp (the exact condition depends on T), the 
right-left equilibration is "bottlenecked" by inelastic col- 
lisions that involve cold electrons near the bottom of the 
electron spectrum. Hence pd acquires the activation fac- 
tor exp(— cf/T) in the low-T limit. The power-law be- 
havior of pd with T that follows from the conventional 
approaches^— is only recovered if the drag rate l/r^ 
[Eq. (|1.6I) ] resulting from the drift ansatz in these ap- 
proaches is smaller than the equilibration rate, propor- 
tional to exp(— cf/T) in the case of soft collisions. At low 
T « tf , this can only occur if the distance a between the 
wires is exponentially large in the parameter ep/T. This 
answers the question formulated between Eqs. (|1.4[) and 
(|1.5[) : the orthodox theory is only justified when drag is 
exponentially weak in £f/T. Conversely, for fixed a, the 
range of applicability of the orthodox theory is limited 
to temperatures which are only "logarithmically smaller" 
than tF- 

On a more detailed note, our results show a nontrivial 
interplay between triple and pair collisions. The activa- 
tion behavior pu oc exp(— e^/T) is determined by triple 
collisions with one electron scattered near the bottom of 
the conduction band and two electrons scattered near the 
Fermi level. If the intrawire triple collisions are less ef- 
fective in the right-left equilibration than interwire pair 
collisions between two cold electrons, there exists a range 
of T in which /?d acquires one more activation factor and 
behaves as exp(— 2ep/T), crossing over into the regime 
dominated by three-electron scattering as T decreases. 
However, in any case, pu is exponentially suppressed at 
low T . One more conclusion that comes from the com- 
parison of the effect that pair and triple collisions have 
on pd is that the orthodox theory of drag for the case 
of forward scattering is totally unjustifiable if only pair 
collisions are present and hinges entirely on the triple- 
collision rate being sufficiently high. 



Our theory of Coulomb drag is built on the quasiclas- 
sical kinetic equation approach. Although we will fo- 
cus most of our attention on scattering with momentum 
transfer much smaller than kp, this approach allows us 
to easily incorporate backscattering near the Fermi level 
as well. Throughout the paper, however, we assume that 
T is still higher than the crossover temperature below 
which interlocked charge-density waves^ r 44 ' 47 induced 
by backscattering are formed. In this paper, we special- 
ize to the case of ballistic quantum wires (no disorder) 
and spinless electrons. 

The paper is organized as follows. Section|n]is devoted 
to Coulomb drag due to pair collisions. In Sec. Ill Al we 
introduce the kinetic equation for a double wire and ob- 
tain the high-frequency drag resistivity. In Sec. IIIB| we 
formulate and solve a model which contains interwire for- 
ward scattering but explicitly forbids backscattering — to 
show that there is no dc drag in one dimension with- 
out backscattering. In Sec. IIIQ we obtain the dc drag 
resistivity induced by pair collisions with small momen- 
tum transfer and demonstrate its activation behavior for 
low temperatures. In Sec. Ill Dt we discuss drag resulting 
from direct backscattering at the Fermi level. Section ITTT1 
deals with Coulomb drag in the presence of both pair 
and triple collisions and emphasizes the important role 
of the latter. In Sec. llHAl we write down the kinetic 
equation that describes triple collisions in a double wire. 
In Sec. IIII Bl we discuss singularities that arise in the 
calculation of the three-particle scattering probabilities. 
In Sec. IIII Ci we describe soft triple collisions within the 
Fokker-Planck approach. In Sec. IIII D| we compare var- 
ious channels of three-particle scattering in the double- 
wire system and identify those that are most important 
for Coulomb drag. In Sec. IIII El we consider the effect 
of triple collisions on Coulomb drag induced by soft pair 
collisions and show that three-particle scattering dramat- 
ically enhances drag at low temperature. Our results are 
summarized in Sec. IIVI Some of the technical details are 
placed in the appendices. 



II. COULOMB DRAG IN ONE DIMENSION: 
PAIR COLLISIONS 

A. Kinetic equation for pair collisions 

Our point of departure is the kinetic equation for pair 
collisions in a system of two spatially homogeneous quan- 
tum wires. In one dimension and for the quadratic dis- 
persion relation, this type of scattering does not affect the 
distribution function if both electrons are in the same 
wire — but does lead to a relaxation of the distribution 
function if electrons arc in different wires. Throughout 
the paper we neglect tunneling between wires, so that 
the exchange processes for electrons from different wires 
are absent. We thus have: 

djAh) - eE a d k jAh) = st ff {/} , (2-1) 



5 



where the symbol a = 1,2 distinguishes wires 1 and 2, St^l/} for the case of pair collisions is given [at the lowest 
E a is the electric field in wire a, and the collision integral (second) order in interaction] for a — 1 by 



Sti 2) {/} = (2^) 2 | ^-j^-J^- - k^Sih + k 2 - fa, fa,) 5(6i + £2 - ei, - 6 2 

x {/i(M/2(fa')[l - /i(fei)][l - /a(*a)] - /i(fci)/ 2 (fc 2 )[l - /i(fa')][l - / 2 (fa')]} , (2.2) 

where ei = fc 2 /2m, etc., and ^(g) is the Fourier component of the interaction potential of electrons in different wires 
with the momentum transfer q [given, e.g., by Vi 2 (q) from Appendix |A~1 throughout Sec. [IT] we omit the subscript of 

V12]. The superscript M = 2 in Stj^j/} in Eq. (|2.2|) is meant to indicate that this is a contribution to the collision 
integral of two-particle collisions (M-particle scattering with M > 2 will be discussed in Sec. IIII|) . For St2{/} in 
Eq. (|2.1[) , the wire indices of f a (k) in Eq. (|2.2p should be transposed (momenta 1 f> 2, 1' f> 2')- The product of the 
delta-functions in Eq. (|2.2[) reduces in the case of quadratic dispersion to 

ffi 

S(h + k 2 - k v - k 2 ,)5{e 1 + e 2 - c v - e a ») = j-r <5(fa - fa,)5(fa - fa) , (2.3) 

\ky - fa| 

which gives 

St< 2) {/} = mj^ ]Vi ^lZt]f {A( fc i')/ 2 (fci)[l ~ /i(fa)][l ~ / 2 (fci0] - /i(*i)/ 2 (fci0[l ~ h(ky)][l - / 2 (fa)]} • 



Below, we will focus on the linear response under the 
assumption that the wires are identical; in particular, 
that their chemical potentials and temperature are the 
same. It is then convenient to represent the solution of 
Eq. (|2.ip in terms of the functions g a {k) as 



U{k) = h + g*{k)Td e h 



(2.5) 



where /x = [1 + e( e_CF )/ T ] -1 is the thermal distribution 
function. Linearizing Eq. (I2.1[) in g CT , we then obtain (in 
the to representation): 



-iwgi(k) - 



n\ eE2k + ( 2 )r l 
-iug 2 (k) - = -st\ '{g} 



(2.6) 



where 



|fc' - k\ 

x [ gi (^+3 2 (fa^<7i(fc)-<? 2 (fc')] (2-7) 



and 



C(fa 



cosh[(e - e F )/2T] 



(2.8) 



By introducing the functions <7±(fc) = [gi(k) ± g 2 (k)]/2, 
we thus have: 



9+{k) = 



e(E l + E 2 )k 1 
2mT — ia; + 



(2.9) 



-i^-(k)- e -^^=s^{g}, (2.10) 



where 

^■f/ic^wn-^)]. 

(2.11) 

One simple result that follows immediately from 
Eq. (|2.10[) gives the drag conductivity 021 (defined as 
°2i = 32/ E\ under the condition that E 2 = 0) in the 
high-frequency limit. Iterating Eq. (|2.10[) in the limit of 
large u once yields Rea2i — —e 2 n/mw 2 T ! §' for wrg° >• 1, 
where 

J_ 2 2 

AmT J 2vr S V ; 

is the electron density in one wire and 

(213) 

The sign 00 is intended to emphasize that the scattering 
rate (|2.13p describes high-frequency drag.-- If the main 
contribution to 1/rg 2 comes from momentum transfers 
with \k' - k\ ~ T/v F , then l/r£° cx T 2 [cf. Eq. ([Eg])], 
while if it comes from backscattering with \k' — k\ ~ 2fa, 
then l/rg° oc T [cf. Eq. (fL~2"1) ]. Naively, one might think- 
in the spirit of the Drude theory or, for that matter, the 
memory-function formalism with the memory function 
expanded to second order in interaction — that 1 / de- 
termines drag also at lj — > 0, with the dc drag resistivity 
Pd oc 1/t^°. As will be seen below, this assumption is 
correct for drag induced by backscattering in the close 
vicinity of the Fermi level; however, it is totally wrong 
for the case of forward scattering. 
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B. Absence of friction from forward scattering 

As mentioned in Sec. HI dc Coulomb drag is absent 
[Eq. (|1.4p ] in the absence of interwire backscattering. To 
see this, consider a model in which the interaction ma- 
trix element does not connect electron states with op- 
posite chirality. In this model, backscattering processes 
both near the Fermi level (momentum transfer about 
2kp) and near the bottom of the spectrum (momentum 
transfer much smaller than kp) are forbidden by con- 
struction. It is important that the model excludes the 
latter possibility as well, because the backscattering pro- 
cesses with small-momentum transfer, while being ex- 
ponentially suppressed for T <C £f, still can lead to a 
"leakage of current" between the subsystems of right- 
and left-moving electrons. For definiteness, let us substi- 
tute for Vk'-k in Eq. (I2.11[) a function of k and k' that is 
proportional to the (9-function of the product kk' , which 
explicitly forbids backscattering^ 



V(k' - fc) -> V(k' - k)e(kk') 



(2.14) 



The model (|2.14p is compatible with those used for study- 
ing Coulomb drag due to forward scattering in Refs. |56| - 
|58| . In the Luttinger-liquid formalism, generalized to the 
finite-curvature case, it corresponds to retaining only the 
9ax tyP e 01 interaction^ 

Within the model (|2.14p . the equation for the distri- 
bution function g_ (k > 0) of right-moving electrons can 
be written in a closed form: 



iujg-{k) - 



e{E l - E 2 )k 
2mT 



9-(k)} , (2.15) 



while the distribution function of left-moving electrons 
cj-(k < 0) is related to g-(k > 0) by the condition 
.?-(— k) = —g(k) which follows from the fact that the 
source term in Eq. (|2.10[) is odd in k. Importantly, the 
collision integral in Eq. (j2.15[) is nullified if g~(k > 0) 
does not depend on k [i.e., g~{k) = const(fc)sgn(fc)]. The 
solution of Eq. (|2.15[) can therefore be represented as a 
sum of two terms, 



flL-(Jfc) = ho + h(k) , 



(2.16) 



where the zero-mode term ho does not depend on k and 
has a pole at oj = 0, 



ho 



A{oj) e(E 1 -E 2 )v F 



-iui + 



2T 



(2.17) 



with a residue proportional to a yet unknown constant 
A(0). The equation for h(k) then reads 



iujh(k) 
lo 2^ 



A(u) 



kp 



e(Ei - E 2 )vf 



2T 



m 



2^(fc'-fc)| 2 



C(k 



\k> - k\ 



[h(k') - h(k)] . (2.18) 



Multiplying Eq. (12.181) by ( 2 (k) and integrating both 
sides over k, we eliminate the collision integral to ob- 
tain the connection between A(lo) and h(k) in a form 
that does not contain the collision kernel explicitly: 



A(u) = 



1 du 
ttv p dn 



•dk 



e(E 1 -E 2 )J Q 2tt 



C 2 h 



(2.19) 



where d/i/dn = 2T '/ J^°(dk/2Tr)( 2 is the inverse com- 
pressibility. For T -C ep, 



1_ dji ^ ir 2 ( T 
itvf dn 8 \cf 



(2.20) 



The difference of the compressibility from 1 / wvp at finite 
T will be important for the calculation of the singular (at 
lo — >• 0) part of CT21. 

From Eq. (|2.19[) . the closed equation for h is written 

as 



j[h(k) -h]- 



e(E 1 - E 2 ) 



2mT 



(k-k) 



w c2(fc0 iv (fc ' 

2ir 



kt 



where 



h=± d Ji 
2T dn 



k' - k\ 



°° rlk 

2tt S 



[h(k') - h(k)] , (2.21) 



- m dpi 
7r dn 



(2.22) 



The solution of Eq. (12.21)) does not contain, by construc- 
tion, a zero-mode part and is regular at uj — > 0. It 
follows from Eq. (|2.19[) . then, that A(oj — > 0) is finite 
(neither diverging nor vanishing) and determined by the 
first term in the square brackets in Eq. (|2.19[) . namely 
A(0) = (l/nv F )d^/dn. The solution of Eq. (j2T5t can 
thus be represented as a sum of the singular (at u) = 0) 
term and the regular term as follows: 



9-(k) - 



A(u>) 



B(w,k) 



e(Ei - E 2 )vf 
2T 



(2.23) 

where A(ui), the "source renormalization" factor B(uj, k), 
and the "memory function" M(uj, k) are all regular at 
lo -> 0, and M(0, k) > 0£1 This form of y_ , together with 
Eq. (|2.9[) for g + , dictates a very special type of behavior 
of the conductivity and resistivity tensors (in the space 
of the wire indices) as oj — > 0, as is seen below. 

The electric current in wire 1 [sign + in Eqs. (|2.24[) 
and (|2.25p below] and wire 2 (— ) is given by 

p f°° dk 

31,2 = ^ C\k) k [g + (k) ± g_{k) } . (2.24) 



Using the relation (|2.19p between the singular (h ) and 
regular (h) parts of g_, the current can be expressed in 
terms of only the regular part as 
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1 



Jl,2 



'dfc 



4m 2 T -iu + J a 2tt 



C 2 fc[(£i + £ 2 )fc±(£i-£ 2 )fc] ± 



The conductivity matrix resulting from Eq. (|2.25[) reads 



a(u) 



e 2 VF 



Ai + A 2 Ai — A 2 



-ilo + Q \X1-X2 Ax+A 2 ) +C ^\ —j [ 



'dk 

2m J 2tt 



1 -1 



C 2 /c (ft - ft) 



where 



and 



Ai 



C(w) 



7T a/i r 00 ^ 



16mT 2 <9n J 2n 



1 <9/Lt 

27tw^ <9n ' 
B(w,fc) 



B(u,k') 



-iu + M(u,k) -iu + M{u,k') 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



For T <C £f, the coefficients Ai, 2 are given by [cf. 
Eq. d230b ] 



where 



A - 1 * 7 T 
Al " 2 " 48 I ^ 



A, 



1 7T f T . 

2 - 16 U ' ' (2 ' 29) 



i.e., the diagonal elements of the first matrix in Eq. (|2.26[) 
are close to unity in the limit of small T, whereas the 
nondiagonal ones vanish as T 2 . Note that the singular 
behavior of the nondiagonal elements is determined by 
the T dependent corrections to the coefficients Ai ;2 in 
Eq. (|2.27p . The function C(u) in front of the second 
matrix is also proportional to T 2 at T — > 0. Indeed, 
the integrals over k and k' in Eq. (|2.28[) are determined 
[because of the factors ((k) and C(^')] ^ ^ ne c l° se vicinity 
of k — k' = kp, while the integrand contains a product of 
two factors each of which is zero at k = k'. The vanishing 
of C{uj) at T — > can also be seen from the sum rule for 
the conductivity (see, e.g., Ref. \Et\), according to which 



dui „ . x e 2 n ( 1 
— Rec(o,) = — ( Q 1 



(2.30) 



independently of the strength of interaction. Equations 
(|2.26[) and (|2.30p . combined together, tell us that 



^ReCH = i(A 1 -A 2 ) , (2.31) 



which, in view of Eq. (|2.29|) , means the T 2 behavior also 
for the integral characteristic of the regular part of <r(w). 

It is instructive to represent the conductivity matrix 
for T <C (f as 



e 2 VF 



1 — 7? Si 



l — 3?7 



2 

t 2 



-iu; + 5-2 



■f Si 

•C(w)E 2 



(2.32) 



the matrices Ei i2 are given by 



Si 









' 1 







I 


, s 2 = ( 







(2.33) 



(2.34) 



and the infinitesimally small frequency shifts iSi and 162 
in the singular terms proportional to Si and E 2 , respec- 
tively, are denoted differently to emphasize the different 
origin of possible damping in the two terms. Specifi- 
cally, the singular term proportional to Ei comes from 
the symmetric (in the wire indices) function g+ whose 
singularity is associated with total-momentum conserva- 
tion. Hence 81 — in homogeneous wires, independently 
of the type and strength of electron-electron interaction. 
In contrast, the singular term proportional to E 2 stems 
from the zero-mode function ho whose singularity reflects 
particle number conservation within each chirality. That 
is, 82 — in the model of only forward electron-electron 
scattering. The last (nonsingular) term in Eq. (|2.32[) is 
the contribution of both ft and the regular part of fto 
[the last term in Eq. (|2.25[) ]: its damping is related to 
the equilibration between electrons of the same chirality 
in different wires. 

The matrix structure of Eq. ([2.32(1 with 81 = S 2 = dif- 
fers in an essential way from that proposed for the same 
case of forward electron-electron scattering in Ref. H3- 
The crucial difference is that the prefactor of E 2 in 
Eq. ([2.32(1 is singular at u> = 0, i.e., behaves in the limit of 
small lj as l/(— iui + S^) with 82 = 0, whereas in Ref. l57l it 
is proportional to l/{—iui + 2/rg 3 ), where the scattering 
rate 1 /rg 3 , describing high-frequency drag, is given by 
Eq. ([2~T3]) . Inversion of the conductivity matrix in Ref. 1 5 71 
yielded a nonzero dc drag resistivity pn = ir/e Vft^, 
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i.e., 1/td [Eq. fjl -6|) ] equal to 1/t£°, which also agrees 
with the result of Refs. I56ll58l . In contrast, the inverse of 
the matrix &(ui) from Eq. (|2.32j) is 



p(w) 



e 2 v F (l-rj)(l-3r)- 2iuC) 



-l + 2r] + iojC 
1] + iuiC 



rj + iuiC 
-l + 2r) + iwC 



which at uj — > gives 



Rep2i ui 



(2.35) 



(2.36) 



That is, in the model of only forward scattering the 
dc drag resistivity po = — P2i(u = 0) is strictly zero 
[Eq. (|1.4j) ]. Note that the diagonal dissipative resistivity 
Repn also vanishes with decreasing uj as uj 2 , similar to 
Rep2i- The coefficient in front of uj 2 is, in both cases, 
proportional to C(0) which scales as T 2 in the low-T 
limit. 

It is worth mentioning that nonzero, in contrast to the 
solution of the kinetic equation, pu in the model of for- 
ward scattering was obtained in Refs. [56T - l58l in two ways. 
In Ref. 57, po 7^ was found as a direct consequence of 
the conjectured Lorentzian shape of the u> dependence of 
021. On the other hand, in Refs. |56|J58|. the same expres- 
sion for pd was shown to result from the one-dimensional 
version of the orthodox theory^— at u> = 0. In particu- 
lar, in Ref. the orthodox theory was cast in the form 
of the drift ansatz. The relation between the two ap- 
proaches and the step in the solution of the kinetic equa- 
tion at which the drift ansatz fails are further discussed 
in Appendix [Bl 

Coulomb drag in the dc limit would only occur if 5% 7^ 
in Eq. (|232"|) , namely 



p(u = 0) 



Ae 2 v F C{0)S 2 + A 2 



(2.37) 



We thus see that the scattering processes that change the 
chirality of electrons — recall that it is these processes that 
yield (5 2 7^ — are the only source of nonzero dc Coulomb 
drag. As already discussed in Sec. U one can distinguish 
two main types of backscattering: in the vicinity of the 
Fermi level and at the bottom of the spectrum. The 
contribution of the former to pd is exponentially sup- 
pressed, as exp(— 4fci?a), if the distance a between the 
wires is much larger than the Fermi wavelength. The con- 
tribution of the latter is also exponentially suppressed, 
as exp(— cf/T) [or exp(— 2cf/T), depending on the pa- 
rameters], if T is much smaller than the Fermi energy. 
It follows that the important parameter that controls 
the relative weight of these two types of backscattering 
in pd is the ratio of a and the "thermal length" vf/T. 
For a ^ vf/T, backscattering with momentum transfer 
much smaller than Uf is expected to provide the main 



contribution to pd- This type of backscattering is dis- 
cussed in Sec. Ill Cl below. 

C. Coulomb drag due to soft pair collisions: 
Fokker-Planck description 

Let us consider the limit in which the characteristic 
momentum transfer in Eq. (|2.1ip is much smaller than 
T/vf- For concreteness, we can think of the interaction 
potential given by Eq. (jA.ll) and simplify Eq. (|2.1ip in 
the limit T/v F > \k' - k\ ~ 1/a [sec Eq. l(A~3j)]. In this 
limit, the collision integral (I2.11[) can be written in the 
differential form: 



in< I d( 2 ()(j 1 
2 \~dk~dk ' 2 
mc 1 d ( 4 dg- 
~4~(?dkV° ~dk 



nC 2 



d 2 9- 
dk 2 



where 



dq 
2^ 



q\\v(q)\ 2 



(2.38) 



(2.39) 



In the limit Ta/vF ^> 1, the scattering rate 1/r-g 3 
[Eq. (|2.13l) ]. which describes high-frequency drag, and 
c are related to each other as follows: 



24 



(2.40) 



Being rewritten in terms of the function 

/- = (fi - h)/2 - g-Td e f T = -9-C74 (2.41) 

(i.e., going back from the "g-functions" to the distribu- 
tion functions / CT ), Eq. ()2 .38[) can be cast in the form of 
the Fokker-Planck equation: 



ejEi - E 2 )( 2 k 
8mT 



dJ^ 2 
dk 



(2.42) 



where the current in momentum space J^(k) [related to 
f-(k) and df-(k)/dk locally — at one point k] is given by 



J^ = -D d J- + f_™ 
dk dk 



(2.43) 



with the fc-dependent diffusion coefficient in momentum 
spa& 



D(k) = ^C 2 (k). 



(2.44) 



The solution of Eq. (|2.42[) in the dc limit can be found 
exactly. At w = 0, J^(k) is obtained by integrating 
Eq. (|2.42l) [assuming that lm\ lJ _j.o(c i - 1 /-) = 0, which will 
be confirmed by the solution]: 
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e(gi - E 2 ) 
8mT 



dpp( 2 ( P ) = -e(E 1 -E 2 ) 



k 2 - k 2 



(2.45) 



The boundary condition used in Eq. (|2.45l) is J (2) (fc ±oo) = 0. Substituting Eq. (|2.45|) in Eq. (12.43j) yields a 
first-order equation for f-(k): 



df_ 

dk 



mT 



tanh 



k 2 - k 2 F 
AmT 



e(£?i-£^) 1 
mc C 2 



fc 2 - A: 2 
1 - taDh, -i^ E 



(2.46) 



which should be solved for the boundary condition f—{k = 0) = 0. The solution reads: 



= e(gi - S 2 ) 2 

TOC 

mc 



dp 



1 — tanh 



p 2 -fc 



4mT 



C 4 (P) 



fc / p 2 -k 

dp exp — 



4mT y C 3 (p) ' 



(2.47) 
(2.48) 



Using the parameter T/e F <C 1, Eq. (|2.48l) can be sim- 
plified to 



/- 



16mc 



(7rmT) 1/2 e 2e -/ T C 2 <i> 



'mT 



(2.49) 

where &(x) = (2/-y/7r) J Q exp(— £ 2 )<it is the error func- 
tion. Equation (|2.49|) is the asymptotically exact expres- 
sion valid for not too large energies e < e* (more accu- 
rately, for f* - e > T), where 



e* = 3e F 



2 T 



(2.50) 



For larger energies e — e* 3> T, it follows from Eq. ()2.48[) 
that /_ falls off as a power law<£2 



/- 



e(gi - g 2 ) r 
2c ft 



(2.51) 



Specifically, for all energies e — e F ^> T, /_ is given by 
the sum of two contributions to the integral (|2.48l) coming 
from \q\ of order (mT) 1 / 2 and from \q\ ~ (2TO6) 1 / 2 : 



f- 



e{E 1 - E 2 ) ( ttT 



1/2 



exp 



2mT 



sgn(fc) 



2(mT/7r) 



1/2 



(2.52) 



The ranges of applicability of Eqs. (|2.49[) and (|2.52l) over- 
lap. For all energies T <C e < e* (which includes mo- 
menta around the peaks of /_ at k = ±k F ), the shape of 
/_ as a function of ft is givenSi simply by £ 2 : 

/- - (^T) 1 / 2 e 2 -/T C 2 sgn(fc) . (2 . 53) 
lomc 

The electric current 

3 -_ = i ,- A ) = -±/|y_, (2.54) 



calculated by integrating Eq. (12 .53[) . reads 
e 2 {E 1 ~E 2 )n 



2m~f 



where 



1/2 

7 = 2c f ) < 



2tF \ P ~2e F /T 
TTT 3 



(2.55) 



(2.56) 



i.e., the difference of the dc conductivities <J\x — <J2i — 
2j-/ (Ei — E 2 ) ~ e 2 v F /ir'-f. Taking into account that the 
sum (Tn + 0"2i = e 2 v F /ir(—iuj + 0) [as it follows from 
Eq. (|2.9[) ] and inverting the conductivity matrix, we ob- 
tain that the dc resistivity matrix 



p(u = 0) = ,00^2 



(2.57) 



where the matrix £ 2 is given by Eq. (|2.34l) . i.e., pn 
— P21 in the dc limit. For the dc drag resistivity pd 
— P2i(w = 0) we have 



PD 



7T7 

2e 2 v F 



ire 
e 2 v F 



2e F 
^T3 



1/2 



-2e F /T 



(2.58) 



Equation (12.58)) is in agreement with the conclusion of 
Sec. IIIBI that backscattering is the only source of dc 
Coulomb drag (and should be contrasted with the result 
of Refs. I56I458I where nonzero po was obtained from for- 
ward scattering in the absence of any backscattering). 
The scattering rate 1/td [Eq. (11.61) ]. which describes 
drag in the dc limit, is seen to be given by 7/2. This 
means that 1/td is much smaller, for T < cf, than the 
scattering rate l/rg° [Eqs. (|2TT3l) . ([2~4TH) ] describing high- 
frequency drag: 



'p 

TD 



24 [ - 

7T 



1/2 



3/2 



-2t F /T 



(2.59) 
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Importantly, the ratio (|2 . 59[) does not depend on the 
strength of interaction, with both scattering rates being 
quadratic in V12. 

It may be instructive to discuss the origin of the T 
dependence in Eq. (|2.58j) in more detail. The factor 
cxp(— 2ep/T) means that the relaxation of /_ at the 
Fermi level in the dc limit is only due to the diffusion 
in energy space which leads to the cooling of an electron 
in wire 1 — starting from the Fermi surface down to the 
very bottom at k = — due to the heating of electrons in 
wire 2, followed by backscattering of the electron at the 
bottom and its acceleration in the opposite direction, ac- 
companied by the cooling of electrons in wire 2. This dif- 
fusion cycle, which amounts to effective backscattering at 
the Fermi level, is bottlenecked by electron-electron scat- 
tering at k — (requiring two holes, one in each of the 
wires, near the bottom) — hence the factor exp(— 2ep/T). 

The above picture is also substantiated by the obtained 
k dependence of the distribution function. Counting the 
exponential factors in Eq. (|2.49[) . we observe that df-/dk 
at k — is proportional to exp(e^/T) and /_ grows with 
increasing |fe| until it reaches maximum at |fc| ~ kp, at 
which point it is proportional to exp(2ef/T). This be- 
havior of /_ means that the characteristic relaxation rate 
for /_ at given k is small as exp(— 2c_f/T) for all mo- 
menta both near the Fermi level and below it down to 
the bottom of the spectrum, i.e., the exponential factor 
in the relaxation rate does not depend on k. Moreover, 
Eq. (|2.53[) says that the relaxation rate at oj = 0, in- 
cluding the pre-exponential factor, is accurately approxi- 
mated in the vicinity of the Fermi level for \k — kp | <C kp 
by a k- independent constant. This constant is precisely 
7 given by Eq. (|2.56[) . For \k — Uf\ "C kF, the r.h.s. of 
the kinetic equation (|2.42l) at u = can thus be written 
as —7/-, with the relaxation rate being independent of 
k and determined by the slowest scattering processes at 
the very bottom. 

It is worth noting that the condition T 3> vp/a, 
which is necessary for the justification of the Fokker- 
Planck description [Eq. (|2.38[) ] in the whole range of mo- 
menta I A: I < kp, is not necessary for scattering within 
the same chiral branch to preserve its diffusive charac- 
ter in energy space for cold particles with |fc| <C mTa. 
If T <C vf/cl, forward scattering with |fc| <C mTa is 
still described by Eq. (|2.38p . In contrast, in the range 
mTa <C |fc| < kp it is modified in an essential way 
by strong asymmetry between hopping up and hopping 
down along the energy axis. Specifically, for a particle 
in this range of fc, the characteristic probability of gain- 
ing energy in an elementary hop is much larger than the 
probability of losing energy. The asymmetry factor de- 
pends on how fast |l^ 2 (<z)| falls off with increasing \q\ for 
\q\a ^> 1 compared to the growth of the thermal fac- 
tor exp[(fc + q) 2 — k 2 ]/2mT. For the case of V(q) from 
Eq. (|A.3|) . the asymmetry factor for Ta/vF <C 1 is mainly 
given 71 by exp[(fc|, — k 2 )/2mT] > 1. The right-left re- 
laxation via multiple scattering with small momentum 
transfer is hindered by the asymmetric hopping. How- 



ever, for T <C vf/cl, direct backscattering with momen- 
tum transfer 2kp becomes important, as will be discussed 
in Sec. HlDl 



D. 



Coulomb drag due to backscattering at the 
Fermi level 



In Sec. Ill CI we have calculated the contribution to pd 
[Eq. (12.581) ] that comes for T ^ vp/a from pair collisions 
with momentum transfer much smaller than kp. Friction 
from this type of scattering has been shown to be medi- 
ated by effective backscattering at the Fermi level, where 
"effective" means that backscattering is a result of the 
diffusion in energy space through the bottom of the spec- 
trum. Let us now turn to the contribution to Pd from 
backscattering at the Fermi level with the momentum 
2kp transferred in one transition. For T -C vf/cl, this 
can be calculated from Eq. (|2.10[) straightforwardly by 
removing the ratio \ V(k' — k)\ 2 /\k' — k\ from under the in- 
tegral sign in Eq. (12.111) and substituting \V(2kp)\ 2 /2kp 
for it. The collision integral (12. lip reduces then to the 
out-scattering term ("relaxation time approximation"), 
which gives 



9- 



e(E 1 - E 2 )k 



with 



2mT(- 



7b = MJr 



7b) 



(2.60) 



(2.61) 



where /3\, — V{2kp)/2 r KVF is the dimensionless ampli- 
tude of backscattering. The simple Lorentzian for the 
lj dependence of g- in Eq. (12.60)) means that the high- 
frequency and dc drag rates for the case of backscatter- 
ing in the vicinity of the Fermi level coincide: 1/rrj = 
I/t^ = 7b/2, in stark contrast to drag induced by scat- 
tering with small momentum transfer [Eq. (|2.59[) ]. The 
result for pd reads^S 



PD 



7T7b 

2e 2 v F 



4tt 2 



T 



VF 



(2.62) 



By comparing the contributions to the dc drag resistiv- 
ity from backscattering at the Fermi level [Eq. (|2.62l) ] and 
effective backscattering due to soft collisions [Eq. (|2.58[) ]. 



pp [Eg. gig)] 

p D [Eq. (EES)] 



V(2ki 



V(l/a) 



Ta 

VF 



1/2 



(2.63) 

one can see that the latter mechanism of backscattering 
provides more friction that the former for— T 3> vf/c 
That is, despite the contribution of soft collisions being 
strongly suppressed compared to Eq. (11.31) . for kpa ^> 1 
there is still a wide range of temperature, vpja <C T <C 
€p, in which soft collisions in a degenerate electron gas 
are more effective than direct backscattering. It is worth 
noting that the main contributions to po only come from 
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backscattering at the Fermi level and from backscattering 
at the very bottom of the spectrum, while backscattering 
at intermediate energies plays no role. Indeed, the expo- 
nential factor exp[— 4|fc|a — 2{ep — e)/T] that describes 
(for e < ep) direct backscattering with momentum trans- 
fer 2|fc| is maximized either at e — €f for T < vp /4a or at 
e = for larger T with a sharp (for kpa ^> 1) crossover 
of width in T of the order of 1/ma 2 . If T <C 1/ma 2 , 
the main contribution to pd comes from direct backscat- 
tering at the Fermi level independently of the parameter 
kpa. 

III. COULOMB DRAG IN ONE DIMENSION: 
TRIPLE COLLISIONS 

In Sees. Ill Bl and lll Cl we have shown that the processes 
of thermal equilibration between electrons with different 
(right-left) chirality are absolutely necessary for the drag 
effect. Further, in Sec. Ill Cl we have demonstrated that 
interwire pair collisions are capable of establishing the 
right-left equilibration and that the resulting drag re- 
sistivity is proportional to exp(— 2ep/T). On the other 
hand, in a single wire, where pair collisions in the ballis- 
tic case do not change the electron distribution function, 
energy relaxation has been known to be associated with 
triple collisions, see Refs. [l! and l74l - [79l for the cases of a 
nondegenerate and degenerate electron gas, respectively. 
In the degenerate case, energy-relaxation processes that 
do not require changing the number of electrons with 
the same chiralit y 75 ' 77 ' 78 are much faster than those that 
do i 74 i 76 ' 79 Specifically, while the former are character- 



ized by scattering rates that are power-law functions of 
the characteristic energy scales, the right-left equilibra- 
tion rate due to triple collisions is proportiona l 74 i 76 ' 79 to 
exp(-ej?/T). 

Despite being exponentially suppressed for T C ep, 
the right-left relaxation is qualitatively important be- 
cause triple collisions can change electric current (in a 
finite-size ballistic wire connected to the leads) only if 
they change the difference between the numbers of right- 
and left-moving electrons! 74 ' 76 ! 79 ' 80 In particular, it is 
the right-left equilibration that determines interaction- 
induced corrections to the conductance and to the ther- 
mopower: in short wires, whose length is smaller than the 
right-left equilibration length, the correction a 74 ' 76 ' 79 are 
proportional to exp(— e^/T) (for a similar consideration 
in the bosonic formulation, see Refs. I8lll82l ). The right- 
left relaxation plays also a role in the transport properties 
of inhomogeneous wires ^ 

Below, we study the contribution of triple collisions 
to the drag resistivity. As already mentioned in Sec. HI 
triple collisions can strongly enhance drag at low T — this 
is precisely because of the right-left equilibration rate due 
to triple collisions being proportional to exp(— ep/T), in 
contrast to exp(— 2ep/T) in the case of pair collisions in 
the double wire. 

A. Kinetic equation for triple collisions 

The contribution to the collision integral St CT for the 
distribution function f a {k\) in Eq. (|2.ip that comes from 
triple collisions reads (for a = 1 in St CT ): 



I 

Stf } {/} - J2 aa , 1 - ± Y^T^ 23V2 ,y W ^' 2 '' 2 > 3 W £ i + £ 2 + e 3 - e v - e v - C3 ,) 

X {/l(M/a(M/<r'(M[l - AMI - fM][l - fa'ih)] 

- h{ki)U{k 2 )f a ,{kz)[l-h(k v )][l-f a {kr)][l-h\ky)]}^ (3.1) 

where the sign ' in means that the summation over momenta goes over (quantum-mechanically) distinguishable 
initial and final states [for ease of presentation, the integration in Eq. (|2.2I) is changed in Eq. (13. ip to the summation 
over quantized momenta: below, L is the size of either wire 1 or wire 2] . The tilde over the summation sign in Eq. (|3.1j) 
denotes one more constraint on the momentum summation — this will be discussed in Sec. IIIIBI 
For spinless electrons, the three-particle state |1, 2, 3) is written as 

|1,2,3)„ = D (fci > fc 2> /5 3 ) , (3.2) 
11,2,3),, = J^e^D^fca) , (3.3) 
|l,2,3) c = -i- e ^3 Dc(A . ljfc2 ) ; ( 3.4) 

where we distinguish three cases: 

(a) : all electrons are in wire 1 , [Eq. (|3.2p ] , 

(b) : electron 1 is in wire 1, electrons 2 and 3 are in wire 2 , [Eq. (|3.3p ] , 

(c) : electrons 1 and 2 are in wire 1, electron 3 is in wire 2 , [Eq. (|3.4p ] , 
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with the normalized three- and two-particle Slater determinants (D Q and T>b,c, respectively) given by Eq. (|C.1I) in 
Appendix [D] The indices a, a' in Eq. (|3.1[i are then grouped as follows: a = a 1 = 1 in case (a) and a = a 1 = 2 in case 
(b). In case (c), identically equal to each other contributions to Sti come from a = l,a' = 2 and a = 2,<r' = 1 [the 
two contributions are weighted with a factor of 1/2 each, as is accounted for by the factor (1 + <5 ff(T /)/2 in Eq. (|3.1|) ]. 
The kernel iiw in Eq. (|3.ip is given in cases (a), (6), (c) by 

«; oAc (l / ,2 / ,3'|l,2,3) = 2n\A% btC (l', 2', 3'|1,2,3)| 2 , (3.5) 

where A" r b (1/, 2', 3'|1, 2, 3) is the irreducible (not factorizable into independent blocks not connected by interaction) 
part of the three-particle amplitude A abc (l', 2', 3'|1, 2, 3). At the lowest (second) order in interaction, the amplitude 
A a>biC (l', 2', 3'|1, 2, 3) is written as 

„/ A { ^ h (V, 2', 3'|4, 5, 6)A£l J4, 5, 6|1, 2, 3) 
A abc 1', 2', 3' 1,2, 3) = V a ' b - cV ' ; / a ' bM ' , (3.6) 

" ^6 ( ei+e2 + e 3 )-(e4 + e 5 + e 6 )+i0 1 ; 

where 

45 c (l, 2, 3|4, 5, 6) = (1, 2, 3|v(zi - x 2 ) + v(x 2 - x 3 ) + v{ Xl - x 3 )\4, 5, 6) a>6 , c . (3.7) 

The potential v(x) in Eq. (|3 . T[) is either vn(i) or v\ 2 (x), depending on whether it couples electrons in the same wire 
or in different wires. 

The sign ' in in Eq. (|3.6p has the same meaning as in Eq. (|3.ip . One can remove the constraint on the allowed 
momenta in Eq. (I3.6[) by introducing the following factors: 

A a = - S T {.••}, A b = -V {...}, A f = -V {...}, (3.8) 

a 6^456 l J ' 2 / -^456 l J ' 2 / -^456 l 1 ' V ; 

where {. . .} is the same fraction as in Eq. (|3.6j) . The factor 1/6=1/3! in A Q comes from the permutations over all 
intermediate states (4,5,6). The factors 1/2 in A b and A c come from the permutations over states (5,6) and over 
states (4,5), respectively. Similarly, the restriction on the summation over momenta in Eq. (I3.1[) can be lifted after 
introducing additional factors in Eq. (|3.1[) . Let us denote st^ 3a ' ) ' ( ' 3b ' ) ' ( - 3c ' 1 the contributions to St^ in Eq. (|3.1[) from 
processes (a),(b),(c). Then, 

s ^ 3a) - 2 -G^)E 231 ^K rr i 2 ^ ■■){•■■}. ( 3 - 9 ) 

St^= 2. (^i)£ 23r2 , 3j Arm... ){..-}, (3.10) 
i 3c) =2 7 ri£ 23i , 2 , 3 ,|4 rr | 2 ^ ••){•••}, (3-11) 



St 

where <5(. . .) denotes the delta- function from Eq. (|3.1I) and {. . .} is the sum of the products of the distribution functions 
as given by the expression in the curly brackets in Eq. (|3.ip . The factors 1/2 and 1/6 in Eq. Q3.9P come from the 
summation over states (2,3) and (1', 2', 3'), respectively. Two factors 1/2 in Eq. (13.101) come from the summation over 
states (2,3) and (2', 3'), respectively. The factor 1/2 in Eq. (|3.1ip comes from the summation over states (1', 2'). Note 
that the combinatorial factors in the collision integrals in Eqs. (|3.9p - (l3.1ip are absent in the formalism of Ref. [74], 
where the three-particle scattering rate (in a single wire) was plugged into the collision integral without an y re striction 
on the summation over the initial and final states [cf. the sign ' in Eq. (I3.ip ], The formalism from Ref. |74j was also 
used in writing the collision integral in Refs. 17514791 The irreducible three-particle amplitudes A" T b c are given by 
Eqs. 1C3J-(1IIZ]) in Appendix O 

Similar to Sec. Ill Al it is convenient to rewrite the linearized kinetic equation that includes triple collisions [Eq. Q3.ip ] 
in terms of the functions g<r(fc) [Eq. (|2.5|) ]. By explicitly separating the zero-mode solution g+(k) [given by Eq. (|2.9p 
independently of the number of colliding particles in the collision integral] , the linearized kinetic equation is represented 
in a closed form for the function g^(k): 

n \ e(Ei~E 2 )h ( 2 ) . , (3) r , 
- lujg-ih) — = st^{ 5 } + st y _'{g} , (3.12) 
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(3) (2) 

where the three-particle contribution st_ {g} to the collision integral depends, similar to st_ {#}, only on g~(k). 
Specifically, st_ {g} is written as a sum of three terms associated, respectively, with processes (a),(b),(c): 

stL 3a) { 5 } = ^ E 2 3 r2 , 3 T a(1, ff i , } 1,2,3 ^ (- • l9-(ki>) + g-(kr) + 9-(k 3/ ) - g.ikr) - 9 -(k 2 ) - g.(k 3 ) } , (3.13) 
st - 36) {5} = \ E 231 ^ ^ (1, '5;51 1 ' 2 ' 3) ^(- • •) l9-(kv) 9-(kv) - <7-(M - + 5 _(fc 2 ) + 5 _(fc 3 ) ] , (3.14) 



4^231'2'3' ^(/d) 

stL 3c) { 3 } = i E^a T^ 1 ^) 1 ' 2,3 ^ !- • ■) l9-(ki>) + MM - <7-(fc 3 <) - 9-(kt) - 9 -(k 2 ) + «,_(*&) ] , (3.15) 



with S(. . .) having the same meaning as in Eqs. (|3.9[) - (|3.11[) and 

W atb<c (l', 2', 3'|1, 2, 3) = 2*143 c(l', 2', 3'|1, 2, 3)| 2 CCfcOC^)^)^)^)^^) (3 lg) 

16 

Note that the contributions to st^ a \g} of the differences <?(A; 2 /) — <?(fc 2 ) and g{ky) — g(fc 3 ) are equal in view of the 
symmetry of the kernel W a (l', 2', 3'|1, 2, 3) = W (l', 3', 2'|1, 3, 2), and the same is true for channel (b). In contrast, 
trading momenta (2,2') o (3,3') for given (1, 1') in channel (c) changes the kernel if Vn(q) ^ V\2{q). If one neglects 

the difference between Vu(q) and Vi 2 (g), the contributions to st_ {fir} of <?(fc 2 ') — g(/c 2 ) and g(ky) — ,g(fc 3 ) cancel 
each other. ^ 



B. Divergencies in the three-particle scattering 
rate 

We now turn to the meaning of the tilde over the sum- 
mation sign in Eqs. (i33|) . (j33| - (i33ll . (l3Tl^ - (l3~T5)) . No- 
tice the energy denominators in the amplitudes A™\ c 
[Eqs. (|C.5I) - (|C.7|) ]: being squared in the collision inte- 
gral [Eq. (|3.5[) or (|3.16[) ]. they yield a singularity in the 
kernel of the collision integral of the type 1/A 2 , where 
A, defined for various scattering processes according to 
Eq. (|C.8[) . is the energy transferred in a virtual transition 
to the intermediate state. The singularity is in general 
not integrable in the sense that the numerator does not 
vanish at A = 0. More specifically, it is not integrable 
if electrons possess spin (or pseudospin, as in the case of 
different wires) — for more details, see Appendix |DJ The 
tilde in Eqs. ( ^ .( ^ - (|33Tl) . ([3^ - (i3~15l) is related to 
the proper handling of the 1 /A 2 singularity in the kinetic 
theory, as explained below. 

Let us first recall the relevant aspects of the many- 
particle scattering problem in the vacuum as we know 
them from quantum mechanical scattering theory. There 
are two conceptually important differences between the 
two-particle scattering problem and the M-particle scat- 
tering problem with M > 2. One is related to the defi- 
nition of the scattering matrix for M > 2. In the former 
case, one can unambiguously define the (exact to arbi- 
trary order in the interaction potential) two-particle scat- 
tering operator whose matrix elements modulus squared, 
taken on the mass-shell in the basis of free (with respect 
to the interparticle interaction potential) states, deter- 
mine the scattering cross-section. In the case of M > 2, 
the scattering states may not be definable as asymp- 
totically free — this happens if particles can form bound 
states in the process of scattering ! 84 ' 85 The M-particle T- 



matrix should then include the scattering states in which 
(a part of) interaction does not disappear at infinity and 
remains important for arbitrarily long times after the col- 
lision event (or, conversely, the bound states may exist 
before the collision event and be excited in the process 
of it). A general formalism that accounts for the proper 
boundary conditions in the M-particle scattering prob- 
lem with arbitrary scattering channels is based on the 
Faddeev equations! 84 ' 85 

For the case of a repulsive interaction (assumed in this 
paper), there are no bound states. However, indepen- 
dently of the sign of interaction, there is another essential 
difference between the scattering problems with M = 2 
and M > 2 — which resembles the one mentioned above 
in that it is also related to scattering processes in which 
interaction remains relevant for arbitrarily long times. 
Specifically, M-particle collisions with M > 2 occur not 
necessarily in a compact region in space and time even 
for the case of a short-range interaction potential. For 
example, three-particle scattering (contributing to the ir- 
reducible part of the scattering amplitude) occurs when 
all three particles are simultaneously within the range of 
the interaction, but it also includes processes which con- 
sist of two consecutive scattering events in one of which 
only two particles interact with each other and the other 
event in which one of those particles interacts with the 
third particlej^— The time separating the two events 
may be arbitrarily long. The 1/A 2 singularity 8 ?' 88 is as- 
sociated with this type of scattering, for which the scat- 
tering probability increases not linearly but quadratically 
in time. 

The 1/A 2 growth of the differential three-particle 
cross-section (defined diagrammatically as a sum squared 
of all connected diagrams for the three-particle scattering 
amplitude at given momenta of the incident and outgoing 
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waves at infinity) as A — > is a no-nonsense singularity 
which requires, however, a proper regularization at the 
point A = 0. Clearly, there arises a question about the 
meaning of the cross-section integrated around A = 0. 
In Rcf. l86l a finite density of particles was introduced to 
regularize the divergency of the three-particle T-matrix 
in the collision i nteg ral. In effect, a similar regularization 
was used in Ref. 1871 . where the quantum kinetic equation 
for triple collisions was derived in terms of scattering am- 
plitudes in a "medium" (the gas of interacting particles). 
In a different approachpSi22 it was pointed out that the 
limit A — > and the limit of the distance sent to infinity 
(in the definition of the T-matrix) do not commute. That 
is, the infinitesimal neighborhood of the point A = in 
the differential cross-section requires delicate handling, 
depending on what quantity is calculated. In particular, 
the implications for the intensity of outgoing waves in 
a three-beam experiment in three dimensions were dis- 
cussed in Ref. |88|. Most importantly, by taking the limit 
of an infinitely large distance after the limit A — > 0, the 
integral of the differential cross-section around the sin- 
gularity was demonstrated to be finite As follows 
from the results of Refs. [86ll87l it is the latter order of 
taking the limits that determines the collision integral in 
the kinetic formulation. Specifically, 1/A 2 in the kernel 
of the collision integral should be regularized at A = 
as the real part of a double pole>2&"— 

1 A 2 -e 2 
|AP ^ (A 2 +e 2 ) 2 ' (3-17) 

which yields a finite result for the integral of the differen- 
tial cross-section over a region that includes A = 0, and 
not as the modulus squared 1/|A + ie| 2 , which would 
give a divergent integral^ The tilde in Eqs. (j3. II) . (|3.9[) - 
(f33l^ . ([3~T3| - ([3~T5|) denotes the regularization rule (|3~T7j) . 

The way the singularity at A = is regularized in 
Eq. (|3.17[) has important ramifications for the structure 
of the collision integral expanded in a series in the num- 
ber M > 2 of col liding particles, St CT = E A / St i- M) - As 
follows from Eq. (|3.17j) . a naive extension of the M = 2 
result, assuming that the kernel of St^ f) for given mo- 
menta is proportional to the modulus squared of the cor- 
responding matrix element of the Af-particle T-matrix, 
is incorrect. The thus defined St^. M ' would be divergent 
for M > 2. As pointed out in Refs. [86l|87l for the case 
of M = 3, the expansion over M contains additionally 
counterterms. Specifically, for M — 3: 

Stf =£{...} = £{...}- if > , (3.18) 

where {. . .} is the contribution to Stj^ of a given set 
of momenta as shown in Eq. (|3.ip . ^{. . .} contains 
1/| A + ie\ 2 (modulus squared of the amplitude) and is 

(3) 

thus diverging, and Jj is a counterterm that cancels 
the divergent contribution to Y]{- ■ •} according to the 
rule (|3.17p . The term —I\ can be considered as remov- 



ing from St \ the contribution of two consecutive two- 
particle collisions separated by an infinite time, so that 
in between the three-particle system returns to the mass- 
shell (the independent two-particle collisions are already 

accounted for by the term St-j/: subtracting the coun- 
terterm thus prevents double counting). 9 — The meaning 
of the tilde in Eq. p. 181) is thus that one should not in- 
clude such "real" states in the summation over virtual 
states in the three-particle scattering amplitudes. Ideo- 
logically, the subtraction of the counterterm in Eq. (|3.18[) 
bears similarity to the treatment of triple collisions in 
classical kinetic theory (§17 in Ref. |92| ). 

The necessity to use the regularization (|3.17l) in the 
collision integral has not been part of the discussion in 
the recent wave of interest in three-particle scattering 
in one dimension^! - — In all these works, the kernel of 
the collision integral is written simply as the modulus 
squared of the T-matrix element. This omission is, in 
fact, only justifiable in the case of spinless electrons [in 
our problem, this corresponds to three-particle scattering 
in channel (a), in which all electrons are (pseudo)spin- 
polarized]: the collision integral in this case does not di- 
verge because of a cancellation between the contributions 
of direct and exchange scattering. In terms of the double 
counting discussed above, the absence of divergencies in 
the spinless case can be understood as a consequence of 
the fact that two-particle collisions do not affect the dis- 
tribution function in a single wire. However, in the case 
of electrons with spin (or pseudospin, as is the case for 
two wires in the drag problem) , it is absolutely necessary 
to specify in what way the singularity at A = should be 
treated [Eq. (|3.17l) ]. because the naive representation of 
the kernel as the modulus squared of the three-particle 
T-matrix element leads to divergency^ The technical 
details of how the three-particle scattering rate behaves 
in the vicinity of the point A = in the drag problem 
are further discussed, in the diagrammatic language, in 
Appendix [Dl where we calculate, as an example, the to- 
tal scattering rate for a (pseudo)spinful particle (for a 
similar calculation in the spinless case, see Ref. |94|). 



C. Fokker-Planck description of soft triple 
collisions 

As already noted at the beginning of Sec. IIIIl three- 
particle scattering may substantially enhance drag for the 
case of soft collisions, when the right-left equilibration is 
controlled by a slow diffusion in energy space. We there- 
fore turn now to a description of three-particle scatter- 
ing with small momentum transfer in terms of a Fokker- 
Planck equation, similar to Sec. Ill Cl for two-particle scat- 
tering. Just as in the case of pair collisions, the Fokker- 
Planck approach is justified if T is much larger than the 
characteristic energy transfer. 

The current in momentum space , induced by triple 
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collisions, is related to sti^{<?} in Eq. (|3.12|) by linearized collision integral [Eqs. (|3.13|l - (|3.15|l ]. the con- 

tributions j( 3a )>( 36 )>( 3c ) to of scattering processes 
( 3 ) 4 dJ^(ki) (a), (b),(c) can be exactly rewritten as 

St ~ M = C 2 (fci) dki (3 ' 19) 

[cf. Eq. (|2.42l) for the case of pair collisions]. For the 



jW(k) = - I ^ / dk'[P a (k',k' + q)g-(k')+P a (k',k' + q)] , (3.20) 



k—q 
k 



where 



J (3b \k) = - P- f dk , [P b (k , ,k / + q)g^(k / )-P b (k / ,k , + q)] : (3.21) 

J-oo Z7r Jk-q 

J {3c) (k) = - [°° p- [ dk'[P c (k',k' + q)g-{k') + P c (k',k' + q)] , (3.22) 

J — oo 27T J k—q 



Pa{kuk v ) = ^ x ~ L V W (l', 2', 3'|1, 2, 3)5(. . .) , (3.23) 



Pb,c{k\ : k v ) are defined similarly, with the numerical coefficient 1/12 in Eq. (|3.23p being changed to 1/4 and 1/2 in 
cases (b) and (c), respectively, and 

Pa(ki, k v ) = i x ^^32,3,^(1', 2', 3'|1, 2, 3)<5(. . .)[ 9 -(k 2 ) ~ 9 -{k 2 ,)] , (3.24) 

P b (h,k v ) = \*\l J2 232>3 W b {l\ 2', 3'|1, 2, 3)<5(. . .)[fl-(*a) - S-(fe') ] . (3-25) 

pfh h \ u/n' o' q'i i o que v 9-( fc 2) - g-(fc 3 ) - g_(fc 2 Q +g-(fc 3 Q ,„ 0fi ^ 

Pc{ki,k v ) = - x ^ i 2^232'3' c( - ' ' |l,2,3)d(. ..) . (3.26) 

The integral over fc' in each of Eqs. (|3.20l) - (|3.22l) is taken over an interval whose width is the transferred momentum 
q. At this point, it is important to realize that it would be incorrect to simply expand in q in the integrands of 
Eqs. (|3.20p and (|3.22|) [channels (a) and (c), respectively] in order to obtain the collision integral in the diffusive limit. 
This is because of the exchange processes in the amplitudes (IC.5|) and (|C.7[) that exchange ky with either ky [channels 
(a) and (c)] or k& [channel (c)]. In these processes, the momentum difference that is small in the diffusive limit is 
&2' — k\ or ky — fci, but not ky — k\ the characteristic value of which is much larger than 1/a. The contribution to 
the collision integral of the processes with small ky — k\ [channels (a) and (c)] or ky — ki [channel (a)] is, however, 
the same as that of the processes with small ky — k\. Therefore, the current J^ 3 \k) is obtained in the diffusive limit 
by expanding in q in Eqs. (|3.20l) - (|3.22l) and multiplying the result by a factor of 3 in channel (a) and a factor of 2 in 
channel (c). More specifically, expanding the products P a ,b,c(k\ k' + q)g-(k') in the integrands to first order in k' — k, 
taking P a ,b,c{k' ,k' + q) out from under the integral sign at the point k' — k, and using the property of the kernel 
P a ,bAk,k + q) 

/oo rk 
dq / dk'P a , b , c (k',k' + q) =0 , (3.27) 
-oo J k—q 

which follows from the vanishing of the collision integral in Eqs. (|3.13[) - (|3.15[) at g-(k) = const(fc) and the condition 
that the current in momentum space is zero at \k\ — > oo, we have 

J&(k) * £>( 3 )(fc)%^ ~ C {3) (k) , (3.28) 
ok 



where 



D {3) (k) = l fp-q 2 [3P a (k,k + q) + P b (k,k + q) + 2P c (k,k + q)}, (3.29) 
2 J 2tt 

& 3 \k)= J^q[3P a (k,k + q)-P b (k,k + q)+2P c (k,k + q)} . (3.30) 
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The collision integral St (3) {/} = -dJ^(k)/dk with J^(k) from Eq. (l3~28)) conserves total momentum and total 
energy. 

Note the absence of a drift component [proportional to g-(k) itself, not its derivative dg-(k)/dk] in the contribution 
to J<> 3 \k) that comes from the terms proportional to P a ,b,c{k' ,k' + q). This is a direct consequence of the exact 
condition (|3.27|) . which can be represented in the limit of small momentum transfer (characteristic q — > 0) as 

JdqqPa tbtC {k,k + q) = Jdqq 2 P atb:C {k,k + q) . (3.31) 

It is also worth noting that, while expanding P a ,b,c(k', k 1 + q)g-(k') ink' — k in Eqs. (|3.20p - (|3.22p yields the first term 
in Eq. (|3.28p . expanding P a ,b,c(k', k' + q) would be beyond the accuracy of the diffusive approximation. Indeed, doing 
so would produce the term dD^ 1 (k) / dk in J^(k), where 

D {3) {k) =i I ^Lq 2 [3P a (k,k + q)-P b (k,k + q) + 2P c (k,k + q)} . (3.32) 



The integrands of P a ,bA k , k + l) [ E qs. ([3724]) - ([3726]) ] for 
scattering k — » k + q of electron 1 contain as factors the 
linear combinations of the differences of the distribution 
functions before and after scattering for electrons 2 and 
3. For q — > 0, the conservation of momentum and energy 
gives two solutions for the pair ky — &2 and k^ — ky. 
either /c2' = and ky = k% or /c2' = &3 and ky = 
k2- In the former case, the expansion of the differences 
g(k2>)—g(k2) and g{ky) — g{k^) around the q = solution 
gives higher powers of q compared to q 2 already present in 
Eq. (|3.32[) . In the latter case, the expansion yields^ the 
factor — &2 , whose characteristic value in the integrand 
at q — > is of the order of the typical momentum transfer 
1 1 a. In either case, adding dD^> (k) / dk to J^ 3 \k) for 
the characteristic transferred energy qk/m much smaller 
than T only gives rise to small corrections to the diffusive 
approximation (|3~251) . Similarly, in Eq. (13"73T))) for (k), 
it suffices to expand g{k2>) — 5(^2) and g(ky) — g{kz] in 
the integrands of P a ,b,c(k 7 k + q) to first order in k 2 > — k 2 
and ky — k^, provided that the characteristic transferred 
momenta are much smaller than T/vf- 

In the derivation of the Fokkcr-Planck equation (|3.28p , 
we assumed that the characteristic change in energy of 
the diffusing electron with momentum k in a single scat- 
tering event is much smaller than T. This allowed us 
to treat g~(k), P a ,b,c(k-> k + l), and P Q ,b, c (fc,fc + q) as 
smooth functions of k on the characteristic scale of q. 
One model in which this condition is satisfied for arbi- 
trary k is that of Vn(q) and Vi2(q) falling off sufficiently 
rapidly as \q\ increases beyond the same characteristic 
scale 1/a <C T/vp. Below, we employ this model for 
estimating the relative weight of various scattering pro- 
cesses in Sec. MI Dl and solving the Fokker-Planck equa- 
tion analytically in Sec. MI El Recall, however, that for 
the case of Coulomb interaction, as can be seen from 
Appendix [M the functions V\\{q) and ^12(9) behave es- 
sentially differently with increasing \q\. Namely, Vi2(q) 
falls off exponentially for \q\a ^ 1, whereas V\\{q) falls 
off only logarithmically (for \q\d 3> 1). In the Coulomb 
case, the Fokker-Planck expansion in Eqs. (|3.20p - (|3.22p 
is justified for Ta/vp 2> 1 for scattering processes that 
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involve the interwire interaction, but is not justified for 
channel (a). Importantly, however, the diffusive char- 
acter of the current J^ a \k) remains intact even in the 
case of Coulomb interaction for |fc| <C kp , as will be 
seen in Sec. MI Dl Electron scattering at the bottom of 
the spectrum (which bottlenecks the right-left equilibra- 
tion) can thus be treated in the Coulomb case within the 
Fokker-Planck approach also for channel (a). The gradi- 
ent expansion of the intrawire contribution to the integral 
term C^(k), on the other hand, would not be justified 
in the case of Coulomb interaction. We will return to the 
Coulomb case at the end of Sec. MI El 

In contrast to the differential Fokker-Planck equation 
for the case of two-particle scattering [Eq. (|2.42p ] , (k) 
from Eq. (|3.28p gives an integro-differential equation — 
because the term dC^(k)/dk in St (3) (fc) is an integral 
of the distribution function. With the integral kernel 
from C^ 3 \k), the equation is not exactly soluble, even in 
the diffusive limits To proceed, we make two approxi- 
mations, one of which is parametrically accurate in the 
drag problem for T/ep <^ 1, the other — for a particular 
relation between the strength of the inter- and intrawire 
interaction potentials (the exact condition will be formu- 
lated in Sec. HUD]) . 



D. Identifying relevant scattering processes 

We now simplify the Fokker-Planck equation in the 
limits mentioned in the last paragraph of Sec. MI CI The 
first step is to realize that, similar to the case of pair 
collisions, the right-left equilibration due to triple col- 
lisions is bottlenecked by the slowing down of diffusive 
motion in energy space around the point k — 0. That is, 
when the equilibration rate limits the drag rate 1/td from 
Eq. (|1.6I) . it is sufficient to calculate J (3) (fc) for |fc| kp 
and the momenta k2 and £3 of two other electrons close 
to the Fermi surface. This separation of scales in mo- 
mentum space is justified in the limit T/ep <^ 1. 

The second step is to compare the contribution to the 
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equilibration rate, induced by three-particle scattering, 
of the region in momentum space in which k 2 and k 3 be- 
long to the same chiral branch (k% ~ k 3 ~ ifcp) and the 
contribution of the region in which k 2 and k 3 are on the 
opposite sides of the Fermi surface (k 2 ~ — k 3 ~ ±kp). 
Specifically, let us estimate the contributions to the diffu- 

(3) (3) 

sion coefficient of a hole with k — > 0, D h and D h ' , 

coming from interactions with electrons on the Fermi sur- 
face with the same (++) or opposite (H — ) chirality. The 
corresponding terms in the electron diffusion coefficient 
D^(k — > 0) are smaller by a factor of exp(— ep/T). 

The conservation of momentum and energy for |fc| <C 
|fc 2 |, \k 3 \ gives 



93 - 92 , 9 - -2^2 



for fc 2 



-fe and 



93 



"92 , 



92" 



(3.33) 



(3.34) 



for ki ~ ^3, where g is the transferred momentum for 
electron 1, q 2 = k 2 > — k 2 and q 3 = k 3 > — k 3 . Impor- 
tantly, while in the former case all three transferred mo- 
menta q, q 2 , q 3 are of the same order of magnitude, in the 
latter case \q\ <C 1 92 1 1 1 93 1 - Specifically, the characteris- 
tic value of 1^1 is of the order of the characteristic value 
°f I92I, I93I ~ min{T/«F, l/a} in the former case and is 



smaller by a factor of T/ep in the latter. This means 
that, for Ta/vp ^> 1, the typical length of an elementary 
step in the diffusion process near k = is of order l/a 
for k 2 ~ — k 3 and of order T jaep for k 2 ~ k 3 . 

To estimate the scattering rate near k — 0, calculate 
first the density of final states 



p{q;k, k 2l k 3 ) 



dq 2 
2vr 



dq 3 S(...)6(q + q 2 +q 3 ) , (3.35) 



where S(. . .) is the delta- function in Eq. (I3.23p . for scat- 
tering 2 — > 2' on the surface in momentum space on which 
the conservation of both total energy and total momen- 
tum is satisfied for given q,k,k 2 ,k 3 . For \k\ -C \k 2 \, \k 3 \ 
and k 2 ~ —k 3 ~ ±k F [Eq. (|3.33j) ]. p(q;k,k 2 ,k 3 ) ~ 
I/A'kvf, whereas the characteristic value of p(o; fc, &2, £3) 
for momenta from Eq. (|3.34p (and \k 2 * — k 3 \ ~ T/vf) can 
be seen to be a factor of e,f/T larger. Next, observe that 
the characteristic width of the integration regions in two 
remaining integrals in P a ,b,c{k, k 4- q), over k 2 and k 3 , is 
T/vf for both k 2 ~ — k 3 and &2 — &3- Now, compare the 
characteristic values of the kernel W a ,f,, c (l/) 2', 3'|1, 2, 3) 
for momenta given by Eqs. (|3.33|) and (I3.34|) . To do so 
in an efficient manner, it is convenient to use the fol- 
lowing properties of some of the fractions that appear in 
A£ r 6c (T, 2', 3'|1, 2, 3) from Eqs. fjCty -QCT ty on the three- 
particle mass-shell: 



(fc 3 < 


- h)(k 3 , 
1 


-k 2 ) 


(fcv 


- k 2 )(ky 
1 


-ki) 




- k 2 ){k v 
1 


-*a) 


(kv 


- h)(k v 
1 


-kt) 




- k 2 )(k 2 , 
1 


-h) 


(k 2 > 


- k 3 ){ky 
1 


-fe 3 ) 


(** 


- k 2 ){k 2 , 


-fci) 




-k 3 )(k 3 , 


-fci) 


(k 3 > 


- fci)(fc* 


-fei) 



(3.36) 

In particular, Eqs. (|3 .36|) show that, for the momenta from both Eq. (|3.33|) and Eq. (|3.34j) . two terms in each of the 
sums almost exactly cancel each other for the case of kpa 1 (exponentially suppressed direct backscattering at 
the Fermi level): the characteristic value of the sum is a factor of kpa smaller than the characteristic value of one of 
the terms. The use of the characteristic values of the momenta in the estimates is justified by the regularization rule 
(f3~T7l) . 

Let us denote (A" r b c )dir the amplitudes of direct scattering in Eqs. (|C.5j) - (|C.7|) . These are associated with the terms 
in Eqs. (|C.5j) - (|C.7| inside the square brackets. For \k\ <C \k 2 \ ~ \k 3 \ ~ kp and k 2 ~ —k 3 [Eq. (|3.33|) ] . they simplify 
significantly: 

rn V 11 (q/2)[V 11 (q)-V 11 (g/2)} 

\ A o.) dil -J^ p > K 6 - 6 ') 

(AV) diT ^ % Fll(g/2) ^f^ (g/2) , (3.38) 
(Airr \ „ rn V l2 (q/2)[V 11 {q) + V 12 {q)-2V 11 {q/2)] 

\ A o Jdir - Jl ^2 ■ ^■' iJ ) 

F 

The main simplification is that the only variable on which the amplitudes in Eqs. (|3.37l) - (|3.39l) depend is the momentum 
transfer q (with the characteristic value of q of order l/a, according to the above). For |fc| <C \k 2 1 ~ 1^3 1 ~ kp and 
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&2 ~ k 3 [Eq. (|3.34p ]. the amplitudes are written as 

m Vn(q 3 ) [Vu (93) - V n (q) ] + gaVufaajVii (?) 



(^ r ) di 



ir - L 2 



hp 



, Ailt \ m V&fo) - ViM VM + 93 n'l (53)^12 (q) 

\ A b Jdir - £2 



(A in ) ~ — MggjXMggl - ^12(93) ] + 53^12(93)^11(5) 

V c /dir £2 



93 ^12(93 ) [ Vii(g) - Vi 2 (g)] 
L 2 q 



hp 



V 



(3.40) 
(3.41) 

(3.42) 



where V t ! ra .,(q 3 ) = dV ar r' (q 3 ) / dq 3 . It is worth mention- 
ing that while the product qzV^ a ,(q 3 ) might seem to im- 
ply that I g 3 1 is assumed to be small compared to the 
characteristic scale on which V a(7 / (q 3 ) changes, taking the 
derivative V^ a ,(q 3 ) in Eqs. (|3.40|l - (13.421) . in fact, only as- 
sumes that the transferred momentum that is small in 
this sense is q, which is a much weaker (in the parameter 
T/e F < 1) condition in view of Eq. (j3~3"4"l) . 

Note that there are two essentially different types of 
strong cancellations between various terms in the deriva- 
tion of Eqs. ([313D]) - ([3H5]> , One, controlled by the pa- 
rameter 1/hpa <C 1, is described by Eqs. (|3.36[) . The 
other, controlled by the parameter T/ep <C 1, is related 
to the destructive interference between two terms given, 
in case (a) , by the first and second lines in Eq. (IC5[) and, 
in case (b), by the first and second lines in Eq. (IC.6[) . 
respectively. Importantly, the latter type of cancella- 
tion does not occur in case (c) [Eq. (|C.7[) ] — because of 
the difference between the inter- and intrawire inter- 
action potentials — which gives rise to the large factor 
q 3 /q ~ k 3 /(k 3 — hy) in the term of (A 1 ")^ on the sec- 
ond line of Eq. (|3.42[) . The singularity of the amplitude at 
fc 3 — > &2' is of the type discussed in Sec. IIII Al and should 
be regularized in the kernel of the collision integral ac- 
cording to Eq. p.!7|) . As follows from this regularization 
rule, the contribution of (AJf r ) to the collision integral 
can be estimated by substituting T jvp as a characteris- 
tic value of the difference £3 — h%i. We see, then, that 
the characteristic value of the amplitude in channel (c) 
contains an additional factor of ep /T 3> 1 compared to 
the amplitudes in channels (a) or (b), so that — unless 
V\\{q) and ^12(9) are very close to each other — channel 
(c) gives the main contribution to the collision kernel for 
the case of electrons 2 and 3 having the same chirality. 

Let us now compare the terms in the diffusion coef- 

(3) (3) 

ficient of a hole with k — > 0, D y h ' ++ and D^',, that 
come from interactions with electrons with hi ~ k 3 and 
&2 — —k 3 , respectively. Each of the two terms is a sum 
of the contributions of channels (a), (b), and (c). It is 
instructive to estimate the relative weight of the six con- 
tributions to the total diffusion coefficient by splitting 
each of them into a product (g 2 )i?/2, where (g 2 ) is the 
average of q 2 in the diffusion process whose elementary 



step is momentum transfer q and R is the characteristic 
scattering rate for these elementary steps. The scatter- 
ing rates include the density of states (13.35)) whose char- 
acteristic value was estimated below Eq. p. 351) to be a 
factor of ep/T larger for the case of k 2 ~ h 3 compared to 
the case of hi — —h 3 . The characteristic values of q for 
|fc| <C h F were discussed below Eq. (|3.34l) and the charac- 
teristic values of momentum differences for electrons on 
the Fermi surface — below Eq. (|3.35[) . Piecing everything 
together, we estimate the scattering rates for k 2 — — k 3 
in channels (a),(b),(c) as 



R 



(3a) 
+- 



V F 


\Vn(l/ay 


4 

( 


'T\ 2 


a 









R 



(36), (3c) 



Vp_ 

a 



V u (l/a) 



-1 2 



Vp 



V 12 {l/a) 



Similarly, for hi — k 3 



„(3a) Vp 

K ++ ~ — 
a 



R (36) Vp 
K ++ — — 
a 



V u (l/a) 



R 



(3c) Vp 



Vp 
Vp 

Vi2(l/a) 

Vp 



Vp 



VniT/aep) 

Vp 

Vi2(T/ae F y 

V F 



TV 

e *7 



(3.43) 



(3.44) 



- , (3.45) 



— , (3.46) 

e F, 



V ll {T/aep)-V 12 {T/aep) 



1 2 



Vp 



(3.47) 



In these estimates, we assume that Tcl/vf > E We also 
assume the most common behavior of the intra- and in- 
terwire potentials, namely (as sufficient conditions) that 
\Vu(l/a)\ > \V 12 (l/a)\ and \V u (T/ae F )\ > \V n (l/a)\. 
Here and in the estimates, the arguments of Vn and V\ 2 
are understood as characteristic scales of transferred mo- 
mentum. The estimate for iiS. [Eq. (|3.47j) ] is written 
under the assumption that |Vn(g) — ^i2(9)|/|Vii(g)| > 
T/tp for \q\ ~ T/aep. Note that i ffij differs from all 
other scattering rates in Eqs. (|3.43l) - (|3.47p in that it does 
not contain the small factor (T/ep) 2 . 

-,(3) 



The hole diffusion coefficient D 



which results 

from triple collisions with electrons of opposite chirality, 



hA 
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is thus estimated, by substituting (g 2 ) ~ 1/a 2 , as 
V&(l/a) 



D 



(3) 
h,+- 



(v F a) 



(3.48) 



Provided Vn(l/a) 3> Vi2(l/a), the main contribution to 
D^]. comes from scattering in channel (a) [Eq. 1(3.430 ]. 

(3) 

The diffusion coefficient D h ++1 associated with interac- 
tions with electrons of the same chirality and character- 
ized by ^g 2 ) ~ (T/ae F ) 2 , is determined by two competing 
terms from channels (a) and (c): 



D 



(3a) VAa/a) 



h,++ 



n ( 3c ) , 



(v F a) 3 

V? 2 (l/a) [Vu(T/ae F ) - V 12 {T/ae F )f 
(v F a) 3 



(3.49) 
T 

(3.50) 



Note that, in channel (a), interactions with electrons of 
opposite chirality are much more effective than with elec- 
trons of the same chirality, because D^'v, in Eq. ([3.490 
has two more powers of the small parameter T/e F com- 
pared to D^\_ in Eq. (pOSf . Therefore, the total dif- 
fusion coefficient can be estimated as a sum of only two 
terms, D { ^ + _ [Eq. (EOSl) ] and D ( £f + [Eq. (l3~50l) ]. Now, 

we observe that the term D^ c ^ + is small in comparison 

(3) 

to D h j in two limiting cases: if the interaction be- 
tween the wires is much weaker than inside the wires, or 
if the two interaction potentials are very close to each 
other. That is, in both limits of a large and small dis- 
tance between the wires, interactions of a hole at the bot- 
tom of the spectrum with electrons of the same chirality 
on the Fermi surface can be neglected. Moreover, in the 
crossover regime, when none of the conditions is satisfied, 

(3c) 

the contribution of D h / , to the total diffusion coefficient 

(3) 

is of the same order of magnitude as that of D h , thus 

not leading to any qualitatively new features, either. 

In Sec. IIIIEI below, we therefore focus on the con- 
tribution of three-particle scattering with two electrons 
having opposite chirality on the Fermi surface. In this 



case, the Fokker-Planck equation, upon substitution of 
Eqs. ((3. 370 - 1(3.390 for the scattering amplitudes, is ex- 
actly soluble for the right-left equilibration rate. 



E. 



Interplay of triple intrawire and pair interwire 
collisions 



Recalling the arguments of the very end of Sec. IHI Dl 
we now consider the right-left equilibration due to triple 
collisions within the framework of the Fokker-Planck 
equation (|3.28|) with the scattering amplitudes (13.370 - 
(|3.390 . These amplitudes correspond to the momentum 
configuration in which two electrons of opposite chiral- 
ity are close to the Fermi surface, while the third elec- 
tron is close to the bottom of the spectrum. We treat 
three- and two-particle soft collisions on an equal foot- 
ing by adding to the current in momentum space (|3.28|) . 
induced by three-particle scattering, the component in- 
duced by two-particle scattering [Eq. ((2.431) ]. One im- 
portant consequence of this is that the mechanisms of 
drag and right-left relaxation, rigidly connected to each 
other in the case of pair collisions in Sec. Ill C[ may now 
be disentangled. To describe the new physics that comes 
about from the interplay of triple and pair collisions, the 
most relevant example is that of drag mediated by pair 
collisions only, with triple collisions occurring between 
electrons all of which are from the same wire [channel 
(a) in the above] . Since drag is only possible in the pres- 
ence of the processes of thermal equilibration between 
electrons of opposite chirality (Sec. IIIB0 . drag may (as 
already noted at the beginning of Sec. Mil) be strongly en- 
hanced by intrawire triple collisions. These do not lead 
to any drag effect directly — but do affect drag indirectly 
by providing an additional channel for the thermaliza- 
tion processes which enhance friction induced by inter- 
wire pair collisions. 

In the limit Ta/v F 3> 1, substituting the amplitudes 
([Q70 - (pT3"9"0 in the kernel of Eq. (j3~2l?f . we obtain the 
terms in the diffusion coefficient D^(k) [Eq. ((3.290 ] at 
|fc| <C k F that result from interactions with electrons of 
opposite chirality in channels (a),(b),(c): 



D (3a \k) 

D {3c \k) : 



C 2 (k) / 




512{irv F ) 3 I 




C 2 W / 


'£) 


512(nv F ) 3 \ 




C 2 (fc) 




1024(7TV F ) 3 





dqq 2 V£(g/2)[V n (<?) ~ ^ll(?/2)] 2 , 



j dqq 2 [V n (q/2)V 12 (q) - V? 2 (?/2)] 2 , 



dqq 2 V 2 2 (qm[Vu{q) + V 12 (q) - 2V n (q/2)] 



(3.51) 
(3.52) 
(3.53) 



Without the factor £ 2 (£;)/4, Eqs. (|3.510 - ((3.53p give the diffusion coefficient for a hole with k —> 0, estimated in 
Eq. (|3.480 . Note that the terms corresponding to channels (a) and (b) have equal contributions of the modulus 
squared of the direct scattering amplitude [Eqs. (|3.4O0 . (I3.410 ]. in which \q 2 \, \qa\ <C k F and \k 2 — ks\ ~ 2k F , and of 
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the modulus squared of the exchange amplitude with |&2' — k 3 \, \k 3 ' — &2I <^ &f and |fe — k 3 \ — 2fep. Altogether, 
taking into account the factors of 3, 1, 2 in front of P a ,b.c{k, k + q) in Eq. p.29[) . the exchange processes thus lead to 
multiplication of the contributions of the processes in which all three transferred momenta q, q 2 , 93 are small compared 
to k,F by factors of 6, 2, 2 in channels (a), (b), (c), respectively. In the limit Ta/vp ^> 1, for the relation between the 
transferred momenta from Eq. (|3 . 33[) . the integral term C^(k) [Eq. (|3. 301) ] reduces to a sum of two terms coming 
from channels (a) and (b), while the contribution of channel (c) can be neglected: 



where 



C (3) (fc) ~ AO) (dg-/dk) 



Di(k) = D^ a \k) - D {zb \k) » 



(3.54) 

(3.55) 
(3.56) 



and k T = fdk( 2 = 8irTdn/dfi ~ 8T/v F . In Eq. ([533]) . the function gt_(fc) only enters C (3) (fc) through the fc- 
independent average p.56|) . which greatly simplifies the solution of the Fokker-Planck equation. 
Let us rewrite the Fokker-Planck equation 



e(£i- E 2 )C, 2 {k)k = d_ 
8mT ~ dk 



D t (k) 



dg-{k) 
dk 



C (3 \k) 



in an integral form: 



9-(k) 



o*AG0 



e(A - £ 2 ) 



1 — tanh 



4mT 



4C (3) (p) 



where the total diffusion coefficient 



D t (k) = D^ 2 \k) + D (3 \k) 



(3.57) 



(3.58) 



(3.59) 



describes both pair and triple collisions. The term describing pair collisions in Eq. (|3.59l) is related to the diffusion 
coefficient D(k) from Eq. (pQ4]) by D^(k) = C 2 (k)D(k)/4. Unlike Eq. (pO?)t for the case of pair collisions [D^ -> 0, 
(7(3) __>, in Eq. p.58[) ]. this is not a solution but an integral equation for g-(k). Substituting Eq. (|3.54p in Eq. (I3.58|) 
and integrating Eq. p.58[) with a factor £ 2 (/e)<9/<%, we have an algebraic equation for (dg-/dk), the solution of which 
gives 



(dg-/dk) = j(E, 



1 — tanh 



k' 2 



kp 



dk 2 , . Dj(k) 
k T C{k) D t (k) 



(3.60) 



Using Eq. (|3.60jl in Eq. (|3.54|) and substituting the thus obtained C (3) (fc) back in Eq. (|3.58|) yields the solution of the 
Fokker-Planck equation in terms of two ^-dependent combinations of the diffusion coefficients, D t (k) and D[(k). 

To characterize the relative strength of two- and three particle scattering, we now introduce two constants (inde- 
pendent of k) T> 2 and T> 3 according to D^ 2 \k) = V 2 C, A {k) and £>< 3 )(fc) = V 3 ( 2 (k), so that 



D t (k) = C(k)[V 2 C{k)+V 3 ] 



(3.61) 



The constant V 2 is related to the constant c [Eqs. (|2~3gii . (|2T4"0|) ] by V 2 = mc/16 = 3fc|/4rg°. Note that V 2 does not 
depend on T, whereas T> 3 is proportional to (T/ep) 2 [Eqs. (|3.51|) - (|3.53[) ]. Similarly, we introduce the constant T> x : 



D i (k)=V i ( 2 (k) 



(3.62) 



The shape of the function g_ (k) in Eq. p.58|) and the resulting resistivity pu will now be parametrized by the "diffusion 
constants" T> 2 , T> 3 , V- u and the ratio T/ep <C 1. Note that the relation between T> 2 and T> 3 can be arbitrary, while 
T>\ < T> 3 [and, depending on the relative strength of channels (a) and (b), 2?i can, in general, be of either sign]. In 
the limit of large separation between the wires, the main contribution to both V 3 and T>i comes from channel (a) and 
T>i ~ T> 3 ; in contrast, in the limit of small separation, T>\ <C T> 3 . 

Let us first calculate the average (|3.60l) . The integral in the denominator of Eq. (|3.60l) can be neglected compared 
to unity in the limit T> 3 <D 2 , while in the opposite limit it is close to unity, which makes the denominator small, if 
Vi~T> a : 



dk V-X 2 



k T T> 2 ( 2 + V 3 



A / 2£>2 
V 3 \ 3V 3 



V 2 -^V 3 



(3.63) 
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The integral in the numerator of Eq. (|3.60l) behaves differently depending on the parameter V 3 e eF ^ T /T> 2 for T> 3 <C ©2, 
which gives rise to three different types of behavior for (dg^/dk): 



(dg_/dk)~- A {E x -E 2 )(^) 



x < 



' (1/4Z> 2 ) [ln 1/2 (P 2 /P 3 ) + (tt 1 / 2 ^)^ 

(l/2? 3 )ln 1 / 2 (2? 3 e eWT /^2) , 
I (e F /T) 1/2 /(V 3 -V i + 2V 2 V i /3V 3 ) , 



X>3 < V 2 e-^/ T , 

V 2 e~^' T « V 3 « ©2 , ( 3 - 64 ) 
V 2 < 2? 3 . 



The range of fc that gives the main contribution to the integral in the numerator of Eq. p.60[) and thus determines 
Eq. (|3T64|) in the three regimes is: |fc| < (2mT) 1 / 2 ln 1/2 (V 2 /V 3 ) for the first term and |fc| < (rnT) 1 / 2 for the second 
term in the first line of Eq. p~64|) . \k\ < (2mT) 1 / 2 \n 1/2 (V 3 e eF / T /T> 2 ) in the second line, and \k\ < k F in the third. 
The logarithmic divergency as T> 3 —¥ in the first line of Eq. (|3.64l) only occurs within the diffusive approximation 
and is cut off^ when the characteristic |fc| ~ (mT) 1 / 2 \n 1 ^ 2 (T> 2 /'D 3 ) becomes of the order of mTa. 

We notice from Eq. (|3.64p that (dg-/dk) contains the large factor e F /T to a certain power (different depending 
on the relation between T> 3 and T> 2 ). This means that, for thermally excited electrons in the vicinity of the Fermi 
surface, the derivative dg-(k) / 'dk from Eq. (I3.58|) is mainly given by the integral term proportional to C^(k) — unless 
T>i/T> 3 is small in the parameter T/e F or, if T> i /T) 3 ~ 1, the ratio V 3 /V 2 < (T/e F ) 1 / 2 e" £F/T . That is, unless the 
above conditions are satisfied, two terms in the total current in momentum space, D t (k)dg-(k)/dk and —C^ 3 '(k), 
almost compensate each other near the Fermi surface. The current in real space, however, is much less sensitive to 
the presence of the integral term in the diffusion equation for three-particle scattering, as will be seen below. 

Substituting Eq. (|3. 581) in Eq. (|2.54[) for the electric current j_ [using the relation (I2.54[) between g_ and /_], we 
obtain j_ as a double integral which is reducible, by integration by parts, to a single one. Using further the relation 
Pp 1 = Aj-/(Ei — E 2 ) between j_ and the drag resistivity pu [see below Eq. (|2.56[t ] . we thus have 



1+e -(e-e F )/T 



J 16tt ./ dk V 2 + V 3 cosh 2 [(e - e F )/2T] 



I + e -0-e P )/T + gv 



(dg-/dk) 



(3.65) 



e(Ex-E2) 

where (dg^/dk) is given by Eqs. (|3.60p . p.64[) . From Eqs. (|3.60p . p.65p . pd is represented as p^ 1 — p^\ + p^\, where 



e 2 T , \l + e-^-^/ T ' 



V 2 +V 3 cosh z [(e-e F )/2T] 
is a direct generalization of the result for pair collisions and 

-1 „ 9T87 , dn f (dg^/8k) \ 2 f [dk 



1 



V 2 + V 3 cosli [(e - e F )/2T] 



(3.66) 



(3.67) 



comes from the integral term of the Fokker-Planck equation. The integral in the curly brackets in Eq. (I3.67P is 
discussed above [Eq. (|3.63[) ]. 

I 



In the limit T> 3 — > and T>i — > 0, we reproduce 
Eq. (|2.58p for pn induced by pair collisions [note that the 
integration in Eqs. (!3.65p . (|3.66D is understood^ as lim- 
ited by |fe| <C mTa]. Now, we observe that the integral 
(|3.65p is determined by e < T in a wide range of the ratio 
r D 3 /V 2 (the exact condition is specified below) — because 
of the exponential functions e _ ^ e_£F ^ T in the numerator 
that rapidly decay away from the bottom of the spec- 
trum. It follows that triple collisions become essentially 
important already for T> 3 3> T> 2 e~ eF ^ T , when the second 
term in the denominator of the integrand of Eq. (|3.65p at 
e = becomes much larger than the first one, i.e., when 
the scattering rate for an electron at the bottom of the 
spectrum is strongly enhanced by triple collisions. 

Inspection of Eqs. (|3.65p - (|3.67p shows that pu is given 



by for all V 3 < V 2 e eF ' T (T/e F ) 3/2 , and for V 3 > 
T> 2 e~ eF/T reads: 




V 2 e~ tFlT < V 3 < V 2 e ep/T (T/e F ) 3/2 . (3.68) 

Recall that T> 3 oc T 2 , so that the pre-exponential factor 
in Eq. (|3T68|) scales with T as T 1 / 2 . Notice that both con- 
ditions on the ratio r D 3 /T> 2 in Eq. (|3.68p are very weak for 
T/e F <C 1. The one that limits r D 3 /T> 2 from above comes 
from a comparison of poi and Pd2 : for larger V 3 /T> 2 , 
Py> is mainly given by PD2- To see this, let us write po2 
for T> 2 <C T> 3 , by substituting the last line in Eq. (|3.64p 
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together with Eq. (l3~63)) in Eq. ([3T67| . The result is 



1 X>3 



PD2 



P 3 -Pi 



2P 2 p 



P 2 < P 3 



e 2 fc F e F Pi V 3P 3 

(3.69) 

Equation (|3.69l) shows that pn 2 can only be smaller than 
Pdi in Eq. (I3.68|) if the wires are sufficiently far away from 
each other, so that P3 — Pi <^T> 2 <C P3. In this limit, 
(P 3 /Pi)(P 3 - Pi + 2P 2 Pi/3P 3 ) 2P 2 /3 and friction 
from /?D2 becomes much larger than that from poi if 
V 2 e e "/ T (T/e F ) 3/2 < P 3 . That is, 



2tt P 2 
3e 2 fc F 



P 2 e eF/T {T/e F f /2 < P 3 . (3.70) 



The behavior of po as a function of P 3 /P 2 with P3 
held fixed is illustrated in Fig. [TJ Note a highly nontriv- 
ial point: even if P3 is entirely due to interactions inside 
the wires and thus does not lead to any drag directly, pr> 
shows a plateau in the dependence on P 2 , i.e., the drag 
resistivity in this regime does not depend on the strength 
of interwire interactions [Eq. (|3.68p ]. In particular, this 
means that varying the distance between the wires in 
this regime does not change pd- On the other hand, pd 
in the plateau regime grows with increasing rate of three- 
particle scattering inside the wire, although this type of 
scattering by itself does lead to any friction between elec- 
trons in different wires. Note also that the width of the 
plateau in the dependence on P 3 /P 2 in Fig.[T]is exponen- 
tially large in the parameter e F /T. The inset in Fig. Q] 
illustrates the behavior of pd with increasing P3 for fixed 
P 2 [from Eq. (j2~5g)) to Eq. ([BTtJg]) to Eq. (|37r0J) ]. 

The T dependence of pd is shown schematically in 
Fig. [3J In a wide range of T, the drag resistivity fol- 
lows the Arrhenius plot with the activation energy equal 
to either e F or 2e F . If P3 3> P 2 at T ~ e F (recall that 
P3 scales as T 2 ), there is a range of T below e F within 
which the T-independent result of the orthodox theory 
[Eq. (I3.70[) ] is valid. As T is lowered, pd starts to be- 
have as e" £F/T [Eq. ([J 



If P 3 < P 2 at T ~ e F , 
there is no room for the orthodox theory for T below e F . 
Instead, pd behaves as e~ 2tF l T [Eq. (|2.58[) ] in a range 
of T right below e F , before crossing over into the e _£F / T 
regime. Eventually, drag crosses over into a low-T regime 
in which it is associated with direct backscattering on the 
Fermi surface [Eq. (I2.6I[) ]. 

Another important point to note is that pd in 
Eq. (|3.70l) coincides with the result of the orthodox 
theory^ for Ta/v F > 1, i.e., Eq. (1X70)) can be rep- 
resented as pd — Tr/e 2 v F T^ , where rp is given by 
Eqs. (|2~T3"]) . ([2"^0"|) . see also Eq. (jlj~9|) . The reason for 
this is that Eq. (|3.70[) describes the limit in which three- 
particle scattering is strong enough to produce the right- 
left relaxation rate that is larger than the drag rate 1/td 
[Eq. (|1.6|) ]. In our formalism, pd that results from the 
application of the drift ansatz is thus associated with 
PD2 [Eq. (|3.67p ]. The condition requiring that the equi- 
libration be sufficiently fast severely restricts the range 
of parameters in which the orthodox theory is valid: for 
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exp(-e F /T) 
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exp(£ F /T) 



FIG. 1: Schematic behavior of the drag resistivity po as a 
function of P3/D2 for fixed D3, where P2 and P3 character- 
ize the strength of two- and three-particle scattering, respec- 
tively. Only the exponential factors are shown in the charac- 
teristic scales on the horizontal axis. Increasing the distance 
between the wires for the case of Coulomb interaction leads 
to a similar behavior of po- In the plateau regime, po does 
not depend on the strength of interwire interactions. The de- 
pendence of po on T in three regimes labeled in the figure for 



Ta/v F > 1: (I) T 



-3/2 -2e f /T 



[Eq. (HMD], (II) T 1//2 e" 



e F /T 



[Eq. (pM)) ]. (Ill) const(T) [Eq. ([370)) ]. Inset: po as a func- 
tion of for fixed T>2, illustrates the growth of pr> with 
increasing strength of interactions inside the wires. The char- 
acteristic scales of T>zjT>2 in the inset are the same as in the 
main figure. 



I \ exp(-2£ F /T) 



Pd 




l/e T 



1/T 



FIG. 2: Schematic behavior of the drag resistivity pD as a 
function of 1/T on the log-linear scale for two distances be- 
tween the wires: larger (solid line) and smaller (dashed). 
Regimes I, II, III are labeled similar to Fig. [T] Drag is hin- 
dered by slow thermal equilibration between two electron sub- 
systems with opposite chiralities, which results in the acti- 
vation behavior of po (regime II: po tx e _e - F/ ' T ', regime I: 
po oc e _2EF//T ). In regime II, po does not depend on the dis- 
tance between the wires. In the low-T regime (labeled with 
ll 2k F "), drag is due to direct backscattering on the Fermi sur- 
face. For T < e F , the conventional contribution to po (labeled 
with "orthodox theory" , regime III) is not suppressed only in 
the case of sufficiently strong interactions inside the wires in 
a narrow range of T right below the Fermi energy e F . 
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given T <C £Fj the orthodox theory is only justified if 
the distance a between the wires is exponentially large in 
£f/T, i.e., if drag is exponentially weak in this parameter 
[cf. the condition in Eq. (|3.70p . where T>2 decreases with 
increasing a, whereas T>z in the limit of large a is due to 
triple collisions inside the wires]. 

One more point worth discussing is the difference in the 
characteristic momenta k that give the main contribution 
to the integral in Eq. (|3.65l) in two transport regimes, one 
described by Eqs. (|2.58p . (l3.68p and the other described 
by Eq. (|3.70p . These are momenta at the very bottom of 
the spectrum, \k\ < (mT) 1 / 2 , in the former case and all 
momenta below the Fermi surface, |fc| < hp, in the lat- 
ter. We emphasize, however, that the distribution func- 
tion f-(k) is sharply peaked at the Fermi surface and the 
integral over k in Eq. (|2.54p — in contrast to Eq. (|3.65l) — 
is determined by \k — kp\ ~ T/vp in both cases. What 
is different between the two regimes from the point of 
view of the characteristic values of k is the range of k for 
the scattering processes that give the main contribution 
to the relaxation rate of the nonequilibrium part of the 
distribution function at \k\ ~ kp. In the case of drag lim- 
ited by the slow right-left relaxation in Eqs. (I2.58p . (|3.68l) . 
this range of k is |fe| < (mT) 1 / 2 , as was already discussed 
in a similar context (for the case of pair collisions) in 
Sec. Ill Cl A subtle difference in the shape of the distribu- 
tion function /_(&) = -g_(fc)/4cosh 2 [(e-e F )/2T] in the 
two transport regimes is that g~(k) ~ const(fc)sgn(/c) for 
all \k\ > (mT) 1 / 2 in the case of Eqs. ([Z55]) . ([3^8l (which 
means that the relaxation time approximation with a k- 
independent relaxation time is asymptotically exact in 
this regime) , whereas g- (k) cx k in the case of Eq. (|3.70p 
(which means a homogeneous shift of the Fermi distri- 
bution in momentum space, i.e., the equilibrium state in 
the moving frame). 

The behavior of p^> with varying T [D3 (T) can be rep- 
resented as T>s(ep)(T/ep) 2 ] and strength of inter- and 
intrawire interactions (T>2 parametrizes the strength of 
interwire interactions and decreases as the distance a be- 
tween the wires is increased, whereas 2?3 in the limit of 
large a is determined by interactions inside the wires) is 
conveniently summarized in the following form: 



PD 



2 7T X> 2 

3 e 2 kF cf 



■ 3 f eF \ 3//2 



e -2e F /T f J 
3/2 



1 • 



e F /T n 



III 



(3.71) 

where regimes I, II, III correspond to those in Figs. Q] 
and [H Note that regime III, when present, is al- 
ways separated by regime II from the u 2kp" -regime 
(Fig. [2]) in the limit of kpa ^> 1, since in this limit 
(1/fcpa) ln[2? 3 (e F )/2? 2 ] is small compared to unity. This 
is because T>3(ep)/T>2 ~ /3fk F a 2 , where /3f parametrizes 
the strength of forward scattering, cf. Eq. (|1 .3[) . As a 



result, for kpa ^> 1, the drift-ansatz regime for the case 
of drag dominated by forward scattering can only be re- 
alized if Ta/vF ^> 1, when the orthodox theory yields 
T-independent drag [Eq. Q3.7ip . regime III; Eq. (14) and 
the plateau regime in Fig. 2 in Ref. 56]. In the oppo- 
site limit of kpa <C 1, drag is determined by backward 
scattering on the Fermi surface (regime "2kp" in Fig. [5]) 
for all T < £f considered in this paper. This implies, 
in turn, that there is no room for Eq. (jl.3l) even if the 
thermal equilibration is strong enough to establish the 
drift-ansatz regime. 

We are now in a position to return to the case of 
Coulomb interaction (recall the discussion at the end of 
Sec. MI Cj) . As we see from the calculation for Ta/vp ^> 1 
[where 1 /a was assumed to be a single scale characteriz- 
ing both functions Vn(q) and Vi2(q), beyond which they 
fall off fast enough to neglect momentum transfer with 
\q\a 3> 1], the main contribution to p^ 1 in the whole 
range of V 3 /V 2 < e £F / T (T/e F ) 3 / 2 comes from p~\, with 
the integral term in the current in momentum space pro- 
ducing only a small correction to pr> ■ The resulting drag 
resistivity is determined by scattering of cold electrons 
with \k\ <C kp- The characteristic energy transfer for 
these electrons is much smaller than T even if that for 
electrons on the Fermi surface is of the order of T, as is 
the case for the intrawire three-particle scattering due to 
Coulomb interaction. It follows that the Fokker-Planck 
description of drag for V 3 /V 2 < e eF/T {T/e F ) 3/2 is also 
accurate for the Coulomb case. However, the diffusion co- 
efficient in channel (a), D^ a \k) [Eq. (|3~2^]) ]. should be 
calculated in the Coulomb case [Vii(<?) from Eq. (|A.2[) ] 
without treating the thermal factors in Eq. (|3.16p as 
smoothly changing functions of q compared to the matrix 
elements, in contrast to Eq. (13.511) . Assuming that T> 3 is 
mainly due to three-particle scattering inside the wires, 
we obtain 



X>3 = 



7r(ln2) 2 T 5 
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2 2 
^pVp 



-) In 2 

V F 



S> 1:1721 



for Td/vp ^> 1, while for smaller T the diffusion co- 
efficient in channel (a) acquires four more powers of 
T. In the opposite limit of fast equilibration in the 
frame moving with the drift velocity, i.e., for *C 
e £J? / T (T/ei?) 3 , the three-particle rate drops out from 
the expression for p^, independently of the character of 
three-particle scattering. Therefore, Eq. (|3.70p describes 
the Coulomb case as well. 



IV. SUMMARY 

We have presented a theory of Coulomb drag between 
clean (no disorder) quantum wires based on the kinetic 
equation approach. One conceptually important aspect 
of Coulomb drag that we have highlighted in this paper is 
an inherent link between this phenomenon and the pro- 
cesses of thermal equilibration. We have demonstrated 



24 



that the dc drag resistivity Pd is exactly zero in the ab- 
sence of equilibration between right- and left-moving elec- 
trons. Another way to state this is that forward scatter- 
ing near the Fermi surface with small momentum transfer 
is not sufficient to produce nonzero drag. 

We have given a detailed discussion of the equilibra- 
tion processes in quantum wires. Crucially, in onc- 
dimensional geometry, the right-left equilibration re- 
quires backscattering — either directly in the vicinity of 
the Fermi surface or via diffusion in energy space with 
small energy transfer in one scattering event. The lat- 
ter type of backscattering is favored if the wires are not 
too close to each other. We have shown that the slow 
diffusion in energy space is bottlenecked by scattering of 
cold electrons at the bottom of the spectrum, as a result 
of which py> shows an activation behavior — in contrast 
to the conventional for the drag effect power-law depen- 
dence on the temperature — with the activation energy 
equal to the Fermi energy tp ov2ep, for the cases of three- 
or two-particle scattering, respectively. We have demon- 
strated a nontrivial interplay between the pair and triple 
collisions; in particular, p^> in a wide range of the param- 
eters of the problem does not depend on the strength of 
interwire interactions, while depending strongly on the 
strength of interactions inside the wires. 
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Appendix A: Interwire interaction potential 



function. Note also that Vi2(x) is not necessarily char- 
acterized by a single scale even if v(r) is a single-scaled 
function. 

For dcfinitcncss, let us consider v(r) in the presence of 
a perfectly screening metallic plate located at a distance 
d from the wires (parallel to them). Then 



vn(x) 



1 



1 



(A.l) 



where d is the "radius of the wire" (which is supposed 
to be the smallest spatial scale in the problem) and the 
dielectric constant of the medium in which the wires are 
imbedded is set to be equal to 1. The Fourier-component 
of vn(x) from Eq. (| A. 1 1) is given by 



Vn(q) 



2e 2 



K (qd )-K (q^/d 2 +M 2 



2e 2 ln 



max{|g|,d x }do 



(A.2) 



where -Ko(x) is the Macdonald function. The Fourier- 
component Vi2(q) of the potential Vi2(x) is given by 



Eq. (|A.2[) with the change d 2 . 



dl 



a 2 in the 



argument of the Macdonald functions and shows the fol- 
lowing behavior depending on whether the distance to 
the gate is larger or smaller than the distance between 
the wires: 



Vi2(q) 



2e 2 ln(d/a) , 
2e 2 ln(l/|a|a) , 
4e 2 d 2 /a 2 , 
e 2 (27r/|g|a) 1 / 2 e-l9l a 

\ _ e -2\q\d 2 /a 



\q\ « 1/d < 1/a 
l/d < \q\ < 1/a 
\q\ < 1/a < l/d 

l/a< \q\ . 



(A.3) 



Note the emergence of the characteristic scale \q\ ~ a/d 2 
in the factor in the last line of Eq. (| A.3|) . One can see, 
however, that the characteristic scale of \q\ on which 
Vi2(q) starts to decay exponentially with increasing \q\ 
is the inverse distance between the wires 1/a, indepen- 
dently of the ratio a/d. 



Let us denote as vn(x) and vi2(x) the potentials of 
interaction between electrons residing in one wire and in 
different wires, respectively. They are expressed in terms 
of the potential v(r) created by a charge in the plane 
of the wires (where r is the distance in this plane) as 
v n(x) = v(|x|) and vi2(x) = v(V x 2 + a 2 ), where a is 
the distance between the wires. The potential v(r) is de- 
termined by the polarization properties of the medium 
around the wires; in particular, by the position and di- 
mensions of a nearby metallic gate. In general, the rela- 
tion between the characteristic spatial scales of v(r) and 
Vi2(x) depends in an essential way on the shape of v(r). 
For instance, if v(r) is a monotonically decaying function 
characterized by a single spatial scale d, the characteristic 
radius of vi2(x) is given by d for d 3> a, while for d <C a 
it may be either larger or smaller than d, depending on 
whether v(r) decreases slower or faster than the Gaussian 



Appendix B: Relation between Eqs. (|2.10[) . (|2.11|> 
and the orthodox theory 

The drift ansatz of Ref. HH [see the discussion in Sec. U 
below Eq. (|1.5[) ] is the result of an extension of the or- 
thodox theory of drag22r— to one dimension. As shown 
in Sees. IIIBI and III CI this approach fails totally in one 
dimension for the description of drag due to forward scat- 
tering. In this appendix, we rewrite the kinetic equation 
(|2.10[) in the form that allows one to explicitly identify 
the approximation that is made in the orthodox theory 
but contradicts the solution of the kinetic equation. To 
this end, let us represent the collision integral (|2.1ip in 
terms of the equilibrium polarization operators 



,' dk 



fr(k + q) - f T (k) 
2vr to + iO - [{k + q) 2 ~ k 2 ]/2m 



(B.l) 
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for two wires [cf. Eq. (|1.5[) ]. whose imaginary parts are 
given by 



Imn(w,g) = - — 



*if + i)-*(T-i 



Using the identity 



Mk')[l-fT(k)} = 



f T (k>) - f T (k) 
l-exp[(e'-e)/T] 



(B.2) Eq. (j2TTj) is rewritten as 



(B.3) 



SM9l ~C 2 Wi 2tt |Jfc'-Jfc| sinh 2 [( £ '- e )/2T] iff - (fc) 
Combining Eqs. (|B.2[) . (IB.4p and changing variables to q — k — k' and oj — e — e', we get 



(B.4) 
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m£ 2 (fc) 



dq \q\\V(q)\ 2 [lmU(oj,q)} 



2tt sinli (w/2T) 



[.-/-(A- -'7I-/A (/.-); ,i( W -^ + |-) . (B.5) 
1 m 2m ' 



Substituting Eq. (|B.5[) into Eq. (|2.10p . multiplying the kinetic equation by ek( 2 /4m, and integrating over fc, we have 
the equation for j_ of the form 



iujj. 



{E 1 -E 2 )n_ e f du' f dq q\V (q)\ 2 [lmU(u' , q) 



2m 



Am 2ir 2ir 



sinli (w'/2T) 



muj q 
~~2 



raw q 



(B.6) 



The result for the dc drag resistivity obtained in Ref. \5Q corresponds to the drift-ansatz replacement 



9- 



9- 



muj q 
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(B.7) 



q 2 J \ q 2 ) enT 

in Eq. (|B.6|) [i.e., g~(k) — > kj-/enT} at ui — > 0. If one employs the drift ansatz ()B.7|) for finite uj as well, this leads to 

e 2 (E 1 -E 2 )n 



J- 



2m(-iuj + 2/rg°) 



(B. 



with 



1 [dco [dq q 2 \V(q)\ 2 [lmIL(uj,q)~ 



8nmT I 2tt I 2ir 



sinri (w/2T) 



(B.9) 



[Eq. (|B.9p coincides with Eq. (|2 . 13[) ] . The Lorentzian shape of the ui dispersion for j_ in Eq. (|B.8[) . with the uj 
independent damping rate (|B.9[) . was posited in Ref. l57l In fact, however, as discussed in Sec. Ill Bl the damping rate 
shows a strong dependence on w, vanishing in the dc limit within the model of Refs. I56M58I . 



Appendix C: Three-particle scattering amplitude 



Explicitly, the normalized determinants in Eqs. (|3.2j) - (|3.4|) read 



B a (ki,k 2 ,k 3 ) 
D 6 (fc 2 ,V 



1 



(6L 3 ) 1 / 2 
1 



e ik 1 x 1 ^ik\x 2 e ik!X 3 



a ik 2 x 3 



e ikzxi e ik$X2 ^ik 3 x 3 



(2L 2 )V2 



e ik 2 x 2 gik 2 x 3 
e ik 3 x 2 e k 3 x 3 



, D c (fci,fc 2 ) 
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^ikxxi „ik\x 2 



,ik 2 xi „k 2 x 2 



(C.l) 
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The matrix elements (13.71) are written as 

1 



^l 15 (1, 2, 3|4, 5, 6) = — {[4 3 ,fc 6 ^ 1 +fc 2 ,fc4+fe 5 ^ii(fci ~ h) + S ku kJk 2 +k 3 ,k 5 +k 6 Vii(k 2 - fa) 



L 

+ fa 2 ,k 5 s ki+k 3 M+k6 V ii{k 3 - fee)] - (fc 2 o fc 3 )} 
-(fci «-> fc 2 ) - (fa «-> fc 3 ) , 



(1,2,3|4,5, 6) 



£ [<5fc 3 ,/C6^fcl+fe2,fe4+fc5 ^12(^1 — ^4) + <5fci,/£ 4 '5fc2+fc3,fe5+fc6^1l(fe ~ fe) 

+ hzfahi+ksM+ke ^12(^3 - ^e)] - ( fc 2 «->• fc 3 ) , 
AW(1,2,3|4,5,6) = — [4 3 ,fc 6 (5 fci+fe2,fc4+fc5^ii( fc i - M + Sk u kJk 2 +k 3 M+k 6 Vi2(k2 - fa) 
+ hs,k B ^ki+k a ,ki+k6 - fa)] - (fa <">• fe) 



(C.2) 
(C.3) 
(C.4) 



(the terms in the third line for Aa are understood to exchange momenta in the whole expression within the curly 
brackets, i.e., for 3 "direct" terms there are 15 exchange terms). The irreducible parts of the amplitudes A a ^^ c 
[Eq. (|3.6[) ] are given by 

1 



A„ r (l ,2 ,3 |1,2,3) - -j-2$k 1 + k 2 +k 3 ,k 1 ,+k 2 ,+k 3l 



where 
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(C.5) 



(fc 2 < O ky) 



(C.6) 



(fa, O fc 2 



(C.7) 



(C.8) 



£2+3-3' = (k 2 + k 3 -k 3 ,) 2 /2m. The sign (fc 2 <-> £3) means 
that only the momenta fc 2 and k 3 are transposed (but fc 2 ' 
and k 3 , are not). The amplitudes A^Jl', 2', 3'|l, 2, 3) 
in Eqs. (|C.5|) - (IC.7[) coincide with those derived in Ref.[7i 



(see also Refs.lZSZSl). It is worth mentioning once more, 
however, that while the amplitude of three-particle scat- 
tering is the same in our work and in Refs. I74ll78ll79l . 
the corresponding contributions to the collision integral 
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are not. This is because the combinatorial factors in 
Eqs. p.9[) - (|3.11[) . necessary to prevent double counting 
of the initial and final states in the collision integral, are 
missing in Refs. [BBS 

Note that if it were not for the difference between 
vn(a;) and vi 2 (a;), the "direct" terms (as opposed to the 
exchange terms) in the amplitudes A" r (l', 2', 3'|1, 2, 3) 
and A™(1', 2', 3'|1, 2, 3) [those shown in Eqs. (ICH . (ICTl) 
with the positive sign] would be expressible as se- 
ries resulting from the cyclic permutations (kiky — ¥ 
k^ky — > k^ky — > kikv). Moreover, the direct terms 
would then become the same in ^4 1 b rr (l', 2', 3'|1, 2, 3) and 
A^ T (l',2',3'\l,2,3). In fact, the whole kinetic problem 
for three-particle scattering of spinless electrons in a dou- 
ble wire would then become identical to that for three- 
particle scattering of spinful electrons in a single wire 
with spin-independent interaction. Our drag problem, in 
which generically vn(x) ^ vi2(x) and the structure of 
A^(l',2',3'\l,2,3) and Aj. rr (l', 2', 3'|1, 2, 3) is therefore 
less symmetric, can be viewed as a generalization of the 
spinful problem in a single wire to the case of Ising-type 
anisotropy of the interaction potential in spin space. 



Appendix D: Cancellation of three-particle 
singularities in one dimension 

As discussed in Sec. IIII B[ one of the important differ- 
ences between two- and three-particle scattering is the 
occurrence of nonintegrable singularities in the differen- 
tial cross-section in the three-particle case. These occur 
if the cross-section is written as the modulus squared of 
(the connected part of) the three-particle T-matrix — this 
would be a straightforward extension of the conventional 
formalism for the two-particle case. In fact, the finite col- 
lision integral that describes triple collisions in the kinetic 
equation contains a counterterm [Eq. (|3.18[) ] that cancels 
the contribution of the singularities. The purpose of this 
appendix is to provide technical details that explicitly 
demonstrate the cancellation between the essential sin- 
gularities in the cross-sections of many-particle scatter- 
ing. Specifically, we focus here on the singular behavior 
of three-particle scattering in the case of one-dimensional 
electrons. 

The amplitude of three-particle scattering in 
Eqs. (|C.5I) - (|C.7P shows a pole-type singularity as a 
function of momenta each time the energy A, endowed 
with indices according to the definition in Eq. (|C.8I) . 
transferred in the virtual transition into the intermedi- 
ate state is equal to zero. In the case of scattering of 
type (a), when all colliding electrons are in the same 
wire, the residue of each of the poles can be shown to 
vanish linearly in A — i.e., the singularity is, in fact, 
absent — provided the initial and final momenta of the 
three-particle amplitude conserve total momentum and 
total energy as A varies. Importantly, the regular 
behavior of the amplitude in channel (a) at A = 
results from a compensation of the direct and exchange 



processes in the residue (for a calculation of the total 
scattering rate in the spinless case see Ref. [94]). For 
the amplitude of three-particle scattering that involves 
electrons from different wires, the compensation is 
not complete — because the exchange interaction in the 
absence of tunneling between the wires is only allowed 
within the same wire — and the residue does not vanish 
(a similar situation occurs for spinful electrons in a single 
wire). Thus triple collisions between electrons belonging 
to different wires yield a nonintegrable singularity in the 
modulus squared of the three-particle T-matrix: at sec- 
ond order in the interaction potential for the amplitude, 
the singularity in the differential cross-section is of the 
type 1/A 2 . 

There is one more important aspect of the divergency 
of the thus defined triple-collision rate that is specific 
to one dimension. The divergency does not rely on 
a particular form of the dispersion law; in particular, 
the singularity is present — and remains nonintegrable — 
in the limit 1/m — > 0. The divergent triple-collision 
rate for electrons with a linear dispersion relation in 
one dimension raises the question as to how the kinetic 
equation approach relates to the Dzyaloshinskii-Larkin 
theore m 40 i 97 which says that, at thermal equilibrium, the 
one-dimensional electron system with a linear dispersion 
relation is exactly described in terms of the random-phase 
approximation. This approximation includes pair colli- 
sions only. That is, according to the theorem (and the 
whole bosonization approach^ for that matter), triple 
collisions are "exactly absent" at equilibrium. The con- 
dition of equilibrium is important; however, the diver- 
gency occurs at the level of the structure of the kernel 
of the collision integral, so that, e.g., the out-scattering 
rate from Eqs. (|C.5[) - (|C.7[) diverges in the linear-response 
limit as well, similar to the inverse lifetime of a particle 
due to triple collisions at equilibrium. Below, we resolve 
the apparent conflict between the Dzyaloshinskii-Larkin 
theorem and the divergency in the three-particle scatter- 
ing channel by calculating the scattering rate at order 
Vf 2 "by brute force" diagrammatically for an arbitrary 
dispersion relation 

The singularity at zero A in the matrix elements in 
either channel (b) or (c) is not related to the difference 
between Vn(q) and Vi2(q) (the singularity survives when 
the difference is neglected) but is only due to the "lack" of 
exchange processes in these channels compared to chan- 
nel (a). Since the singularity is entirely associated with 
scattering of electrons belonging to different wires, we 
neglect intrawire interactions throughout Appendix iDl 
Moreover, since our purpose in this appendix is to il- 
lustrate the principle (discussed in Sec. IIII A[) on which 
the cancellation of the 1 / A 2 divergencies is based, we do 
not calculate here the full set of out- and in-scattering 
nonequilibrium self-energies for two- and three-particle 
scattering but focus on the simplest quantity that exem- 
plifies the problem. This is the inverse electron lifetime 
in an equilibrium electron bath, expanded in V±2(q) to 
fourth order. In this calculation, the inverse lifetime will 
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a) 



b) 




FIG. 3: (a) Self-energy of the fourth order in the interwire interaction in channel (c) for direct scattering processes, (b) 
Effective interaction (thick wavy line) expressed in terms of the bare interwire interaction (thin wavy lines). The electron lines 
for different wires are labeled by the upward and downward arrows. 



be seen to be a well-behaved scattering rate that experi- 
ences no infrared divergency from the vicinity of the point 
A = 0. The quantities of interest are thus the inverse life- 
times for an electron with momentum k in channels (b) 
and (c), 



lAv(fc) = -2ImS fc:C (ie„ -> £ fc + iO, k) 



(D.l) 



where S(, iC (ie n , k) are the corresponding electron self- 
energies in the Matsubara representation, at order V-f 2 
and zeroth order in V±i- 

One can separate the contributions of direct (H) and 
exchange (F) processes in the self-energy in Eq. (ID. 1 1) . 



J b,c 



yH 

^b,c 



yF 



(D.2) 



For the case of triple collisions, the H-term comes from 
the sum squared of the terms with sign + in Eq. (|C.6j) 
or (|C.7|) [for channel (b) and (c), respectively] plus the 
sum squared of the terms with sign — , while the F-term 
is given by twice the product of the two sums. In fact, 
the self-energy Sb jC contains also a contribution of pair 
collisions at order V^ 2 , for which one can similarly sepa- 
rate the direct and exchange processes. The role of pair 
collisions will be discussed below Eq. (|D.26|) . Since in 
channels (b) and (c) the H- and F-terms do not compen- 
sate each other, it suffices — for the purpose of describing 



the divergency of the triple-collision rate — to focus on 
one of the terms: below, we write down details of the 
calculation for the if-term only. 

Let us begin with channel (c) — by calculating the scat- 
tering rate for an electron in wire 1 due to interaction 
with two other electrons, one of which is in wire 1 and 
the other is in wire 2. The self-energy of fourth order 
in interaction for the case Vu = is given by the dia- 
gram in Fig. [3^, where the thick wavy line is the effective 
interaction V(ifl m , q), shown in Fig. [5Jd and written as 



V(in m ,q) = \V 12 (q)\ 2 Tj2 



dq' 
2^ 



W)l s 



x A{ifl m ,q\iQ m ,,q')Il(i£l m ,,q') . (D.3) 

Here the polarization operator Tl(ifl m ,q) is the general- 
ization of Eq. (|B.1[) to arbitrary , 



Il(i£L m ,q) = 



dk f T (k + q) - f T (k) 



(D.4) 



and A(iQ m , q\i£l m r, q') is a sum of the four-leg loops in 
Fig. 2]over all nonequivalent insertions of one of the legs: 



f dk 

A(in m ,q\in m ,,q') = -Tj2 / ^ 



1 



+ 



(^n £,k ) (J'^n ^k — q'} (^n i£l m ' £,k — q' ) (^n ^^m— m' £,k — q — q' ) 

l 



{it-n - £,k) 2 {ie n - i&m - £,k-q)(ien - i&m> - Cfe-g') 

1 



£,k-q-q' ) 



(D.5) 



For the case of a linear dispersion relation, A vanishes (apart from the uncertainty at ift m — > + iO and similarly for 
ifl m >) in agreement with the loop cancellation (Dzyaloshinskii-Larkin) theorem.] 40 ' 97 However, as already mentioned 
above, the divergency of the integral of the modulus squared of the three-particle T-matrix elements over A around 
A = survives the limit of the linearized dispersion law. 
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FIG. 4: Sum of the four-leg loops from the effective interaction in Fig. [3b. The legs are labeled with the incoming frequencies. 



Summing over e„ in Eq. (ID.5|) and Q m > in Eq. (|D.3I) . X^(ir2 m ,(?) reads 

va^ q) = |y 12 ( g )P / g / g [M*0 - M« + <rtl / g + ^ _ & 

x {[/T(fc - 5) - /r(fc - </)] + [/t(* -q-q')~ h(k - q)} z 2 + [f T (k) - f T (k - q)] z 3 

dfrjk) c(6 - 6-9O - c(£k>+g> - 60 dfrjk - g) c(6- g - 6-g-gO - c(6'+g' - 60 
96 Ai 96- 9 A 2 

where c(w) = (1/2) coth(w/2T), 

Ai = 6 + 6' - 6-9' - 6'+«' ) ^2 = 6-9 + 6' - 6-9-9' _ 6'+9' i 

and the functions 21,2,3 are given by 



(D.6) 



(D.7) 



z\ = 



+ 



Ai 



dc(oj) 



c(6 - 6-9') ~ c (6'+9' ~ 60 c(6-9 - 6-gO - c(£k>+q> - 6') 



Ai 



iCl m + 6-9 — 6-9' — 6'+9' + 6' 



1 c(6-« - 6-<?0 - c(6 - 6-«0 



Ai 



ifl m + 6-9 — 6 



1 



c(6 - 6-9' ) - c(6'+9' - 6' ) c(6-9 - 6-9' ) - c (i 



k'+q' — 6' 



Z2 



1 

A^ 



dc(oj) 



duj 



+ 



Ai id m + 6-9 — 6-9' — 6'+9' + 6' . 

c(6 - 6-9' ) - c (6'+9' - 6' ) c (6 - 6-9-9' ) _ c (6'+9' - 6' ) 



(D.8) 



1 C(6 - 6-9-9') ~ C (6'+9' ~ 6' 
\ 2 i^m + 6-9 _ 6 



Ai 



-i£i m + 6 - 6-9-9' — 6'+9' + 6' 



1 



c(6-9 _ 6-9-9') _ c (6'+9' - 60 c(6-9 - 6-9') _ c (6'+9' - 6') 



-iQm + 6 ~ 6-9-9' — 6'+9' + 6' 



(D.9) 



Z3 



1 

aT 



i9c(w) 



9uj 

1 



c(6 - 6-9') - c (6'+9' - 6 



Ai 



c(6 - 6-9O - c (6'+9' - 60 c (6 - 6-9-9O _ c (6'+9' - 60 



i^m + 6-9 — 6 

1 

i^ljn ~t~ 6—9 — 9' 6- 



Ai — iQ m + 6 — 6—9-9' — 6'+9' + 6' 

c(6-9 _ 6-9-9O _ c (6'+9' - 60 c(6 - 6-9-9O _ c (6'+9' - 60 



-if2 m + 6 ~ 6-9-9' ~ 6'+9' + 6' 



.(D.10) 



The terms with dfr(k) /<96 = ~C 2 (' C )/47 1 in Eq. (|D.6I) arise from the double poles in the fermionic summation over 
e„. The terms with dc(w)/duj = — 1/4T sinh 2 (w/2T) in Eqs. ()D.8j) - (|D.10|) arise from the double poles in the bosonic 
summation over fl m >. The terms in the second lines in Eqs. (|D.8I) - (|D.10|) are proportional to (if2 m + 6-9 — 6) _1 an d 
produce, together with the same factor in the first line in Eq. (|D.6|) . double poles in £l m . 
The self-energy in Fig. reads 



£f(k„,feO = -T^ J 



V(iVt m ,q) 
2n ie n - i£l m - 61- 



(D.ll) 
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Doing the summation over fl m in the terms of S c that are proportional to 21,2,3 gives 



1 

a7 



9c(w) 



9w 



+ 



"(6 ~ ~ c(6'+g' ~ 6') 

Ai 



J(«en,6i-gl6-g _ 6) 



' g ^ q ^2 ^ +g [-^ e m6i-gl6-g - 6) _ I(itn,£,k 1 -q\£,k-q + 6' _ 6-9' _ 6'+g')] 

c(6-g - 6-9') - c(6 ~ 6-9') dl{ie n ,£ kl - q \v) 



duj 



1 C (6 6-9') c (6-9 6-9O [r/- elf f \ T{; c f If CM 

H 7— 7T 7T\ l-'(. ze n,6i-gl6-g-g' - 6-9') _ 1 ( le n, 61-9 1 6-g ~ 6)J 

Ai(A 2 — Aij 

+ 9 ^ Q } A ^ +g — [I fan, 61 -g I6-9-9' - 6-9') ^ ^(«e«,6i-9l6-9 +6' - 6-9' - 6'+g')] , 



Z2 



(z6 n zXi m ^/ Cl _g)(ir^ m -|- ^k — q £fe) 



1 

A2 



9c(w) 



9w 



+ 



c (6-9 ~ 6-9-9') ~ c (6'+g' ~ 60 



6i-g|6-g - 6) 



+ 



c (6 6-9-9') c (6'+g' 6') t It t \ TC- t It if t t 
^2 m ze «;6i-9l6-9 _ 6) ~ J-^n, 6i-gl6-g-g' + 6'+g' ~6 ~ 6', 

C(6 - 6-9-9') - c(6-g - 6-9-9') ^(^n,6i-gk) 



<9w 



^> — Ifc-g — ?fc 



c (6-9 6-9-9') c (6'+g' 6') r r /. f \ t if c \i 
A 2 (A 2 - Ai) [J(^n,6i-gl6-g-g' - 6-9') ~ 1 ( ie n, 6i-g|6-g - 6)J 

. c (6 — 6-9-9') — c (6'+9' — 6') r r /. t It t \ r/„-, t it t M 

H Ai(A 2 - Ai) [J(^n,6i-gl6-g-g' - 6-9') " J ( le n, 6i-g|6-g - 6)J 

1 C (6 _ 6-9-9' ) _ C (6'+9' — 6' ) r T 1 ■ >. |> l-\ T(- <- It it I- t 

H -7— \ 61 -g I6-9 - 6) 6i-9l6-g-g' + 6'+g' - 6 - 6' 



Z3 



(ie n - iQ m - 6i-g)(^m + 6-9 _ 6) 
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aT 



dc(co] 



duj 



+ 



c(6 - 6-9O — c (6'+g' ~ 6') 



-f(»en,6i-gl6-9 _ 6) 



H — — ^ 9 9 ^ 2 ~ +g ~ [I (ten, 61-9I6-9 _ 6) - J(«en,6i-gl6-g-g' + 6'+g' - 6 - 6')] 



1 1 
AT ~ A^ 



[c(£fc - 6-9-9') - c (6'+9' _ 6' 



dl(ie n ,£ kl - g \u) 



duj 



c(6 6-g') c (6 6-9-9') c If C \ J li c tit CM 
A^A 2 - Ai) [ J (^n7 6i-9l6-g - 6) - -<^en,6i-gl6-g-g' - 6-g')J 

c (6 ~ 6-9-9') ~ c (6'+g' ~ 6') r r /- t If C ^ T(; c f If _L f C C 

7— T ?fci-g|?fc-g ~ 4fcj - J U e n , 4fci -9 ls/s-9-9' + 4fe'+g' ~ ?fe ~ ?fe' 

AiA 2 



where 



I(te n ,€\u) 



c(u)-t(Q 
te n - £ + w 



and = (l/2)tanh(e/2T). 

The imaginary part of the retarded self-energy comes from the functions I: 

7T 

lml(ie n -> £ kl + i0,6i-g|u>) = -- [c(<; fel _ 9 - ^ fel ) - i(6i-g)j ^(61 - 61-9 + w ) • 



(D.12) 



(D.13) 



(D.14) 



(D.15) 



(D.16) 
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Triple collisions are associated with the terms in Eqs. (ID.12[) - (|D.14[) that contain the functions I with six electron 
energies in the denominator [i.e., six electron energies in the delta-function in Eq. (|D.16[) ] . namely I(ie n , £fci-g|£fc-g + 
£k' - €k-q> - Cfe'+g') and I(ie n ,£, kl -q\£k-q-q> + £,k'+q' - ik — 6c')- These are only present in the contributions to 
coming from the functions £1,2,3- More specifically, they are absent in the combination z 2 + z 3 , so that regrouping the 
terms proportional to #1,2,3 in Eq. (|D.6|) as 

[Mk -q)- f T (k - q')] Zl + [f T (k -q-q')- fr{k)]z 2 + [f T (k) - f T (k - q)](z 2 + z 3 ) , (D.17) 

only the first two differences fr(k — q) — fr(k — q 1 ) and fr(k — q — q') — fr(k) describe the rate of triple collisions. 
Thus we obtain the contribution l/r ( ? 3 (fci) to the triple-collision rate in channel (c) from direct scattering: 

^ [%z \vMi)\ 2 \ d 4z IW)I 2 [c(^-q - &J - ] [Mk' + q ') Mk')] 



r "s(ki) J 2tt 1 KHn J 2tt 1 XH n J 2tt 

[ c(&_, - ) - c(&/ +9< -&,)][ Mk -q)-Mk- q') ; 

x <K6vi + Ck-q + 6e' - £k 1 -q - £k-q> - £,k'+q>) 











/ 2vr 






A 2 ) 



+ A~ 2 ^" aTJ [c(Zk-Zk-q-q')-c(Z k , + q,-{; k ,)]lMk-q-q , )-Mk)} 

X 5(£ kl + ik-q-q> + £k'+q' - &CX-9 - tk - 6')} • (D.18) 

Two different delta- functions in Eq. (|D.18[) correspond to two different amplitudes: the one proportional to A^ 1 (A^ 1 — 
A^ 1 ) comes from the process shown in Fig.[SJi, the other — from the process shown in Fig. [5}). The energy denominators 
Ai and A2 for the processes in Figs.JS^, and[SjD are then identified in terms of the energies in Eq. (|C.8|) as follows: 

Fig. 3a: Ai = -A 3 n< , A 2 = A 322 < , (D.19) 
Fig. 3b: A 1 = A 3 ir, A 2 = -A 322 < . (D.20) 

The sum of the contributions of the two processes to Eq. (|D.18[) gives the cross-section proportional to 

2 



1 1 



(D.21) 



A311' A3 22 ' 

[cf. Eq. (|C.7|) ]. This is because the factor 

MZ kl -q - 6=J - *(&!-,)][/(*' + q') - f(k')][f(k -q)- f(k - q')Mt k _ q - ^ k _ q ,) - c(£ k , +q , - &,)] (D.22) 

from the process in Fig. [5^ and its counterpart from the process in Fig. [SJd, 

[c(£ kl -q - & J " t^ kl -q)}[f(k' + q') - f(k')}{f(k -q-q')- /(fc)][c(& - - c(£ k , +q , - &,)] , (D.23) 

become identically equal to each other when written in terms of the energies £1,2,3,1', 2', 3' in Fig. [SJ Specifically, each 
of them is written as 

— [coth(l - 1') + tanh(l')][tanh(2) - tanh(2')][tanh(3) - tanh(3')] [coth(3 - 3') + coth(2 - 2')] , (D.24) 
16 

where coth(l — 1') = coth[(ei — e±t)/2T], etc. Further, for t\ + e 2 + £3 = £1/ + £ 2 < + £3', Eq. (|D.24I) can be reduced to 

cosh(l) , D 

16cosh(2)cosh(3)cosh(l')cosh(2')cosh(3') ' 1 ' ' 

which is recognized as the factor that appears in the collision integral (|3.1[> for the scattering rate 1 — > 1', namely 

/(*a)/(fcs)[l - f{kv)][l - f{k 2 ,)][l - f(k 3 ,)} + f(k v )f(k 2 ,)f(k 3 ,)[l - f(k 2 )][l - f(k 3 )} , (D.26) 

taken at thermal equilibrium. We have thus reproduced the inverse lifetime due to triple collisions diagrammatically. 

The scattering rate in Eq. (|D. 181) diverges because of the factors 1/Af and 1/A|. But, as can be seen from 
Eqs. (ID.12[) - (|D.14[) . the factors 1/Aj and 1/A| are present in the self-energy (jD.ll) not only in the part associated 
with triple collisions [six fermionic energies in the argument of the delta-function in Eq. (ID.16[) ]. but also in the part 
that contains the delta-function of a sum of four fermionic energies and is therefore identified with the contribution 
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FIG. 5: Three-particle scattering amplitudes at second order in the bare interwire interaction (wavy lines). They contribute to 
the first (a) and second (b) terms in the kernel of the triple-collision rate in Eq. (|D.18[) . The electron lines for different wires 
are labeled by the upward and downward arrows. 




FIG. 6: Cuts labeled 2-2 and 1-3 in this particular term in the three-particle self-energy contribute to the total three-particle 
scattering rate and to a reduction of it that comes from two consecutive two-particle scattering events, respectively. Each of 
the contributions is diverging (with opposite signs), their sum is finite. 



of pair collisions at order V^ 2 - These scattering processes are associated with the product of two amplitudes in which 
one is of order V12, the other — of order V^- For example, if one considers the diagram for the three-particle self- 
energy in Fig. [6l the two-particle processes correspond to the cut labeled 1-3 (one amplitude is of the first order 
in interaction, the other — of the third order), in contrast to the cut labeled 2-2 which contributes to l/r ( f f 3 (fci) in 
Eq. (ID. 181) . Specifically, we have for the sum R c ^{ki) of the terms in l/r c (fci) that originate from direct scattering 
in pair collisions and are proportional to cither 1/Af or 1/A|: 



RcMh) = 2n J^-\V 12 (o)\ 2 J^\V 12 ( q ')\ 2 J^[c(^ 



dk 
2^ 



+ 
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1 

Af 
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c(£k -6-?') - c(6s'+g' -&')] 

■ Cfc-g-g') - c (6 - £,k-q-q 



T2W6 



£,k — q — q' 



X <J(£ fcl +£k-q - ffcx-Q 



IfT(k-q)-Mk-q') 
f T (k) - f T (k - q)} 
)][h{k-q-q') 
f T (k) -f T (k- q) 



Mk-q)} 



(D.27) 



which diverges in the same manner as l/r^ 3 (ki). To see how the sum of the two contributions to l/r c behaves, 
compare the factors in front of 1/Af in Eqs. (ID. 181) and (ID.27[) . Importantly, the delta- functions 5(£ki + 6c-q + 
6c' - £fci- g - - £k>+q') in Eq. (|D.18I) and S(^i + Zk-q - Zki-q - Ik) in E q. (|D 27|) become identical at Ai = 0. 
Further, the difference c(£& — — c (£fc'+g' — £fc') in the third line of Eq. (ID.27|) vanishes at Ai = 0, while the 

remaining factor in front of 1/Af in Eq. (|D.27|) and its counterpart in Eq. (|D.18I) exactly cancel each other. A similar 
cancellation occurs with the terms proportional to 1/A|. We thus arrive at the conclusion that, in channel (c), the 
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FIG. 7: Two contributions to the self-energy of the fourth order in the interwire interaction (wavy lines) in channel (b) for 
direct scattering processes. The electron lines for different wires are labeled by the upward and downward arrows. 

divergency in the total triple-collision rate, defined via the integral of the modulus squared of the matrix elements of 
the three-particle T-matrix, is exactly cancelled by the divergency associated with pair collisions. How this happens 
diagrammatically in the above calculation of the inverse lifetime at equilibrium exemplifies the subtraction of the 
countcrterm (|3.18[) in the three-particle collision integral. 

Let us now turn to channel (b), where the self-energy E^ (ie n , ki) is given (at order Vyx and zeroth order in Vu, 
similar to the above) by the sum of two diagrams in Fig. [Jj 

Ef(ie„,fci) = T 2 j^\ V M? /^|V 12 ( ? 0| 2 5(ie„,fcx|^ m , 9 ; ifi mS «0n(^m, 9^(^,9') , (D.28) 

mm' 

where 



Bfe n ,ki\iCl m ,q;iSl m >,q') = 



1 



i^m £fei— q)\^T 
1 



1 



26 n iil m _ 



(D.29) 



As far as the singularities are concerned, the structure of E^ is much simpler than that of E~% because the noninte- 
grable singularities of the type 1/A 2 2 encountered separately for triple and double collisions in Eq. (|D.11[) are absent 
in Eq. (|D.28j) . Specifically, the only^S term in E^(ie n , fci) that contains a factor similar to 1/A 2 2 in channel (c) is 

^ \y^(l)\ 2 I % I W)| 2 / g IMk) - f T (k - q)} [f T (k> + q') - f T (k')} 



X ^2 I C (^i-9 ~ &!-«-«') _ C (Cfc'+ 9 ' - Zk')][Ifen,€ki-q-q>\€k-q + £,k' - 6c - tk>+q>) - I fen, €ki-q\£k-q + 6)] 



(D.30) 



where A3 = + £,k'+q' — q — £,k' ■ However, in contrast to channel (c), the difference of the bosonic 

distribution functions c(^ 1 _ g — £ki-q-q') — c{£,k'+q' — £,k') in Eq. (ID.30[) vanishes at A3 = 0, so that the singularity 
reduces to I/A3 and the integral (taken in the principal value sense for the scattering rate) is finite. 

Note that, at order and zeroth order in Vu, the electron self-energy contains more terms than the diagrams with 
two electron loops, shown in Figs. [3] and [71 and their exchange counterparts (with one loop less). These come from 
further rcnormalizations of the two-particle T-matrix (with one loop and three loops) not included in the exchange 
counterparts of Figs. [3] and [7J in particular, from the diagram with a chain of three loops, which corresponds to the 
effective interaction in the random-phase approximation. Summing up contributions of two-particle cuts of higher- 
order self-energy diagrams will give the modulus squared of the matrix elements of the exact (in V12 , in this particular 
case) two-particle T-matrix. For the precise meaning of an M-particle cut, see Ref. [87l 

As seen from the above calculation, the divergencies in the triple- and pair-collision contributions to the scattering 
rate and their cancellation occur for an arbitrary form of the dispersion relation of colliding particles. In particular, 
this means that the nonintegrable singularity in the differential cross-section for triple collisions that comes from the 
modulus squared of the three-particle T-matrix is present in the Luttinger liquid model (linear dispersion) as well. The 
three-particle singularity is, however, canceled by the other one that comes, in the kinetic-equation formalism, from 
those two-particle collisions in which a given particle participates twice. This resolves the apparent conflict, mentioned 
at the beginning of this appendix, between the divergency of the triple-collision rate and the Dzyloshinskii-Larkin 
theorem for the linear dispersion law. 
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In Ref. l7q . which deals with triple collisions in a short 
(not equilibrated) single wire [channel (a) in our nomen- 
clature], the validity of the Fokker-Planck approximation 
for Coulomb interaction was argued to be controlled by the 
parameters T/ef <C 1 and |fe| <C fci? — since the character- 
istic transferred momentum is then given by T/vf for each 
of three colliding electrons and the characteristic change of 
energy for the diffusing electron is T\k\/kF <C T. However, 



the integral term (k) was missing in the Fokker-Planck 
equation written in Ref. [7f3 [because the nonequlibrium 
corrections to the coefficients A and B in Eqs. (B.5),(B.6) 
in Ref. [76| were not included in the collision integral] . The 
omission of C (3) (fc) may strongly affect the relaxation rates 
and, moreover, violates the conservation of momentum. 
The condition for the Fokker-Planck expansion in the inte- 
gral term C*' 3 ^ (k) is different from those specified in Ref.l76l. 
namely — since the distribution function in C^ 3 '(k) is taken 
in the vicinity of the Fermi level — the characteristic trans- 
ferred momentum for electrons with velocity Vf should be 
smaller than T/vf- This requires that Vn(q) be a sharply 
peaked function on the scale of T/vf, which means that T 
is limited from below, in contrast to Ref. [76|. 
I.E. Dzyaloshinskii and A.I. Larkin, Sov. Phys. JETP 38, 
202 (1974). 

It is worth noting that E^(ie n ,fci) contains other terms 
that have a transferred energy squared in the denominator; 
however, these have the structure of either dl(ie„, /dui 
or dl(ie n ,^\uj) /d^, none of which leads to a divergency in 
ImSf (ie n -> f fcl +i0,ki). 



